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I. 


INTRODUCTION 


In  combustion  testing,  it  has  been  difficult  to  evaluate  the 
mechanisms  of  erosion  because  of  an  inability  to  separate  the  effects  of 
pure  forced  convective  heating  of  the  bore  from  combined  thermochemical 
effects.  Use  of  propellants  containing  carbon  and  oxygen  as  a  means  of 
providing  bore  surface  heating  and  erosion  is  clearly  in  conflict  with  the 
desire  for  separation  of  effects,  but  has  until  now  been  the  only  means  by 
which  heating  conditions  of  sufficient  magnitude  could  be  produced.  With 
the  introduction  of  the  Shock  Tube  Gun  (STG),  described  in  this  report, 
Calspan  has  provided  a  tool  which  overcomes  the  limitations  of  the  combus¬ 
tion  approach  to  erosion  testing,  and  can  permit  examination  of  the 
thermochemical  nature  of  erosion  in  the  near  melting  range. 

Prior  research^  developed  the  measurement  techniques,  for  use  with 
the  STG,  which  permit  accurate  determination  of  material  loss  and  heating. 
Concurrent  formulation  of  a  computer  code  which  translated  ballistic  and 
gross  heating  data  into  bore  surface  temperature  histories  and  melting  rates, 
predicted  the  heating  conditions  required  to  produce  pure  melting  of  steel 
in  the  STG.  Testing  corroborated  the  melting  predictions  of  the  code  and 
thus,  the  primary  objective  of  the  previous  study  was  achieved.  However, 
computer,  STG  test  analysis,  and  60mm  firing  tests2  conducted  at  Calspan  have 
indicated  that  significant  erosion  and  cracking  may  be  caused  by  chemical 
reactions  in  propellants.  This  report  describes  the  results  of  ongoing 
research  to  identify  and  quantify  the  thermochemical ly  induced  erosion  and 
cracking  in  gun  barrels. 

The  primary  objective  of  this  study  was  to  quantify  the  erosion  of 
a  typical  gun  steel,  SAE  4340,  in  a  gun-like  environment  using  gas  mixtures 
that  varied  in  their  oxidizing  potential.  Therefore,  the  amount  of  material 
loss  for  given  test  gas  conditions  was  a  primary  measurement  of  causality  in 
the  role  that  oxygen  plays  in  gun  barrel  erosion  and  cracking.  This  primary 
objective  has  been  largely  achieved. 

A  secondary  objective  was  to  develop  the  capability  of  predicting 
erosion  in  terms  of  measurable  test  quantities.  A  predictive  computer  code 
was  developed  which  simulates  the  thermodynamic  mechanism  and  governing 
relationships  in  the  Shock  Tube  Gun  during  a  test.  Ultimately,  this  code 
will  provide,  in  correlation  with  empirical  findings  of  the  STG,  a  model  of 
propellant  erosivity  and  the  ability  to  minimize  erosion  features  in  future 
propellant  formulations.  Presently,  this  portion  of  the  study  requires 
refinement  and  remains  an  area  for  further  research. 


1.  F.A.  Vassallo,  W.R.  Brown,  "Shock  Tube  Gun  Melting  Erosion  Study,"  USA 
Ballistic  Research  Laboratory,  ARRADCOM,  Contract  Report  ARBRL-CR-00406 , 
January  1979.  (AD#A076219) 

2.  F.A.  Vassallo,  "Study  of  Wear  and  Erosion  in  the  60mm  MC-AAAC  Gun," 

Calspan  Monthly  Reports  to  ARRADCOM,  Nos,  1  through  12,  to  1  October  1978. 


II.  SHOCK  TUBE  GUN  FACILITY 


2.1  STG  Concept 

Cal span  has  developed  a  unique  variable  gun  system,  the  operation  of 
which  is  based  upon  shock  tube  principles  (£hock  Tube  Gun).  In  contrast 
to  measuring  erosion  that  results  from  firing  a  vast  number  of  combustion 
rounds,  the  facility  is  designed  to  be  used  as  a  developmental  tool  for 
propellant  formulation  changes  and  to  assure  an  approach  to  optimal  propellant 
tube  interfaces.  In  achieving  this  goal,  tests  in  the  Shock  Tube  Gun  are 
supported  by  special  erosion  and  thermal  sensors,  metallurgical  analysis, 
and  ballistic  measurements.  Briefly,  the  STG  generates  a  high  pressure, 
high-temperature  test  gas  by  an  adiabatic  compression  process.  As  shown 
in  Figure  1,  the  facility  consists  of  a  driver  chamber,  driven  tube,  a 
flying  piston,  a  gas  collection  chamber  and  an  instrumented  gun  tube  con¬ 
taining  a  projectile.  Compression  of  the  test  gas  (the  counterpart  to  the 
propellant  gas  in  an  actual  gun)  within  the  driven  tube  is  accomplished 
by  motion  of  the  flying  piston  contained  within  this  tube.  Motion  of 
the  piston  is  the  result  of  force  applied  by  the  pressure  of  the  driver  gas. 

As  the  piston  approaches  the  collection  chamber,  the  test  gas  pressure 
and  temperature  rapidly  increase  in  a  time  history  representative  of  actual 
gun  firings.  A  projectile  contained  at  the  start  of  the  gun  tube  is  acted 
upon  by  this  collected  gas.  Suitable  means  are  provided  in  the  facility 
to  allow  regulation  of  shot  start  pressure  on  the  projectile.  Once 
released,  the  projectile  is  accelerated  along  the  tube  in  a  ballistic  cycle 
dependent  upon  selected  input  factors.  Data  regarding  pressure,  velocity, 
heating,  and  erosion  are  collected  through  measurements  in  the  instrumented 
test  barrel. 

With  suitable  variation  in  parameters,  factors  affecting  erosion 
such  as  pressure  history,  propellant  gas  velocity,  gas  temperature,  gas 
composition,  tube  composition,  and  propellant  additives  may  be  investigated. 
The  compression  ratio,  driven  gas  composition  and  its  initial  conditions 
essentially  govern  the  peak  temperature  and  pressure;  piston  motion, 
influenced  by  its  mass,  and  projectile  movement  effectively  govern  the  pres¬ 
sure  pulse  duration.  Variation  of  piston  mass,  initial  conditions,  and 
compression  ratio  then  permit  independent  change  in  peak  temperature  and  pres¬ 
sure  as  well  as  time.  Obviously,  effects  due  to  variation  of  driven  gas  com¬ 
position  may  be  tested  under  controlled  interior  ballistic  conditions.  This 
is  a  most  powerful  experimental  mode  of  operation  of  the  facility,  and 
the  principal  mode  used  in  this  study. 


2 . 2  Design  and  Construction 

The  Shock  Tube  Gun  is  designed  to  adequately  represent  predicted 
ballistic  conditions  within  guns  as  large  as  the  eight-inch  howitzer, 
requiring  peak  test  chamber  pressures  up  to  300  MPa  and  projectile  velocities 
up  to  800  m/sec.  Table  1  lists  the  present  structural  makeup  of  the  Shock 
Tube  Gun.  These  values  were  established  through  the  use  of  a  preliminary 
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mathematical  model  of  piston  action  based  upon  adiabatic  compression.  Through 
exercise  of  the  model,  approximate  size  requirements  were  established  with 
final  selection  of  sizes  dictated  by  available  engineering  materials.  The 
photographs  of  Figure  1  illustrate  the  resulting  Shock  Tube  Gun  assembly 
as  well  as  its  individual  components.  As  shown,  the  projectile  launch  com¬ 
ponents  consist  of  the  driven  tube,  chamber  and  30mm  smooth  bore  barrel. 

These  are  supported  on  a  carriage  which  is  free  to  move  on  tracks  in  the 
direction  of  piston  motion,  during  the  extreme  impulse  loads  imposed  by  the 
unbalanced  chamber  pressure  during  test  gas  compression.  This  maintains  the 
integrity  of  the  supporting  base  structure.  The  floating  mount  system 
requires  use  of  a  pneumatic  brake  on  the  driven  tube  for  safety. 

The  projectile  capture  components  consist  of  a  telescoping  tube 
coupled  to  the  barrel,  a  blast  chamber,  and  a  sand  filled  tube  to  decelerate 
the  projectile  (impact  zone).  The  blast  chamber,  the  purpose  of  which  is 
to  reduce  the  noise  and  pressure  levels  at  projectile  exit  of  the  tube,  also 
contains  an  internal  provision  for  measurement  of  projectile  velocity  using 
velocity  screens.  The  telescoping  tube  allows  motion  of  the  carriage  inde¬ 
pendent  of  the  blast  chamber. 

Also  shown  in  Figure  1  are  views  of  the  chamber  and  toggle  restraint 
system  needed  to  contain  the  high  chamber  pressures  and  associated  axial 
loads.  Chamber  pressures  are  sensed  using  piezoelectric  pressure  transducers. 
The  entrance  region  of  the  launch  tube  can  accommodate  pressure,  heat  flux, 
and  erosion  sensing  devices.  The  launch  tube  itself  is  a  30mm  smooth  bore 
barrel,  4.57  meters  long. 

Another  essential  component  of  the  facility  is  the  brake  area  at  the 
upstream  end  of  the  driven  tube.  This  brake  limits  the  maximum  permissible 
axial  load  on  the  driven  tube  during  the  rapid  piston  deceleration  period. 
Without  it,  loads  could  exceed  the  axial  strength  capability  of  the  tube. 

In  essence,  load  is  limited  by  slippage  in  the  brake  at  a  preselected  load 
below  the  failure  strength  of  the  tube.  The  brake  is  air  actuated  and  thus 
may  be  adjusted  according  to  the  amount  of  slippage  desired  up  to  the  strength 
limit  of  the  tube.  The  brake  is  presently  used  at  only  about  one-half  capacity 
without  excessive  slippage.  Hence,  much  greater  maximum  chamber  pressures 
than  presently  produced  can  be  accommodated. 

The  piston  which  is  used  to  compress  test  gas  in  the  driven  tube  is 
made  from  4340  steel  and  weighs  68  kg.  inc'  ing  the  release  catch  at  its  rear 
end.  Gas  seal  is  obtained  using  "T"  rings  at  the  front  and  rear  of  the  piston. 
Three  brass  wear  rings  are  used  to  prevent  metal-to-metal  contact  between  piston 
and  tube.  A  buffer  projection  on  the  face  of  the  piston  and  a  complementary 
close-fitting  port  at  the  end  of  the  driven  tube  entering  the  chamber  preclude 
direct  impact  of  the  piston  with  the  end  of  the  driven  tube  in  the  event 
that  driven  gas  is  exhausted  too  rapidly,  as  by  a  failure  in  the  chamber. 

The  buffer  projection  is  necessitated  by  the  presence  of  the  chamber  volume 
at  the  end  of  the  driven  tube  which  can  permit  piston  contact  at  the  driven 
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tube  end  with  insufficient  chamber  pressure  development.  Metal -to-metal 
impact,  if  allowed  to  take  place,  could  cause  damage  to  the  piston  and/or 
tube  face.  The  buffer  prevents  this  occurrence. 

For  any  particular  combination  of  driven  (test)  gases  and  desired 
ballistic  result.  Shock  Tube  Gun  settings,  namely:  driver  pressure,  piston 
and  projectile  weights,  and  driven  gas  initial  pressure  are  determined  using 
a  computer  code  which  models  the  dynamic  conditions  within  the  gun.  The 
computer  code  was  validated  in  early  developmental  testing  of  the  STG 
whereby  piston  position,  chamber  pressure,  and  projectile  velocity  were 
measured  and  compared  with  those  predicted.  Using  the  computer  code,  STG 
input  factors  can  be  chosen  to  result  in  preselected  test  conditions  of 
pressure,  temperature,  gas  velocity  and  gas  activity  at  the  entrance  of  the 
launch  tube.  These  represent  parameters  which  are  important  in  tube  heating 
and  erosion. 


Table  1 .  Shock  Tube  Gun  Characteristics 


Configuration  Data: 


Driven  Tube  I.D. 

0.191m 

(7.5  in.) 

Driven  Tube  Length 

24.6m 

(970  in.) 

Piston  Area 

0.0285m2 

(44.179  sq.  in.) 

Piston  Mass 

Up  to  91  kg. 

(200  lb.) 

Projectile  Diameter 

30mm 

(1.181  in.) 

Projectile  Area 

706mm^ 

(1.095  sq.  in.) 

Projectile  Mass 

Up  to  .91  kg 

(2  lb.) 

Driver  Volume 

0 .885m'5 

(54,000  cu.  in.) 

Chamber  Volume 

2140mm3 

(130.8  cu.  in.) 

Pressure  -  at  release  of  projectile 

Tube  Length 

Variable 

4.57m 

(180  in.) 

III.  TEST  PREPARATIONS 


3.1  Specimen  Selection  and  Installation 

The  chief  objective  of  this  work  is  to  determine  propellant  gas 
conditions  which  lead  to  possible  oxidation  and/or  cracking  of  the  bore  sur¬ 
face,  resulting  in  material  loss  and  fatigue.  To  do  this  effectively,  test 
bore  specimens  amenable  to  analysis  had  to  be  uniform  in  shape  and  composition, 
to  eliminate  these  variables  in  the  testing,  and  representative  of  a  conven¬ 
tional  gun  steel.  4340  steel,  which  has  been  increasingly  used  in  gun  barrel 
applications  since  the  early  1950*  sf/ was  chosen  as  the  test  bore  material. 

The  physical  properties  of  4340  steel  are  listed  in  Table  2  .  Tests  conducted 
with  4340  in  the  previous  program*  failed  to  create  a  melting  condition. 
Therefore,  the  test  sample  was  altered  by  decreasing  its  bore  diameter  from 
20mm  to  12.7mm  (0.5  inches)  .  This  flow  channel  reduction  creates  a  sonic 
flow  condition  which  elevates  the  surface  heat  flux.  The  test  sample  design 
is  shown  in  Figure  2.  Also  shown  in  the  figure  are  the  ports  for  in-wall 
thermocouples  which  are  used  to  determine  the  integrated  heat  input. 

The  erosion  test  section  within  the  chamber  of  the  Shock  Tube  Gun 
is  as  shown  in  Figure  3.  Each  test  specimen  is  fabricated  in  the  form  of  a 
cylindrical  nozzle  and  placed  downstream  from  the  entrance  throat.  The  outer 
diameter  of  the  specimens  was  selected  to  be  31.75mm  and  each  was  38.1mm  long. 
Samples  of  this  size  fit  on  the  specimen  stage  of  Calspan's  Etec  Autoscan 
Scanning  Electron  Microscope  (SEM)  such  that  "before  and  after"  examination 
of  the  bore  surface  can  be  conducted  without  recourse  to  replicas. 

Each  test  specimen  was  characterized  with  regard  to  mass  and  bore 
diameter  before  testing.  Specimens  were  weighed  using  an  analytical  balance 
so  that  mass  changes  of  a  fraction  of  a  milligram  could  be  readily  determined. 
The  bore  diameter  was  measured  to  within  10_4mm  at  four  specific  axial 
locations  as  measured  from  the  leading  edge  of  the  sample:  6.35,  18.4,  27.3, 
37.8mm  (0.25,  0.725,  1.075,  1.49  inches). 


3.2  Heat  Transfer  Instrumentation 


A  primary  measurement  of  the  study  is  the  amount  of  bore  heating 
associated  with  each  test.  For  this  measurement,  in-depth  thermocouples  were 
installed  in  selected  samples.  These  in-depth  thermocouples  were  placed  at 
distances  approximately 0.5mm  from  the  bore  surface.  The  method  of  installa¬ 
tion  is  as  shown  in  Figure  4.  Each  of  these  thermocouples  independently  may 
be  used  to  determine  net  heating  to  the  bore.  Total  heat  input  is  calculated 
from  the  in-wall  thermocouples  output  by  use  of  methods  developed  and 
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Figure  3  SHOCK  TUBE  GUN  TEST  SECTION 


Figure  4  IN-WALL  THERMOCOUPLE  INSTALLATION 


reported  by  Calspan.^  Briefly,  conversion  of  in-wall  thermocouple  outputs 
(millivolts  vs.  time)  to  total  net  heat  input  per  unit  area  is  made  by  use 
of  the  relation 

Q(t)  ■  AT(t)  Akcpt  (1) 


where  Q(t)  is  the  net  bore  heat  input 

AT(t)  is  the  indicated  change  in  in-wall  teaperature  as  a  function 
of  time 

k.  is  the  thermal  conductivity 

cp  is  the  heat  capacity  per  unit  volume 

t  is  the  time  after  start  of  heating. 

Data  reduction  procedure  is  simply  to  apply  Equation  (1)  at  successive 
time  intervals  (e.g.,  0.05,  0. 1, .0.15  sec,  etc.)  resulting  in  a  plot 
of  Q(t)  vs.  t.  The  curve  thus  produced  will  be  asymptotic  to  the  desired 
heat  input. 


3.3  Selection  of  Test  Conditions 


The  program  objective  was  to  experimentally  characterize  the  role 
of  oxidation  in  gun  barrel  erosion.  Oxidation  reactions  in  guns  are  produced 
by  carbon  dioxide  and  water  vapor.  These  reactions  are  suppressed  by  the 
presence  of  carbon  monoxide.  It  is  believed  that  the  neutral  point  occurs 
when  the  ratio  C0/C02  equals  3.  The  combustion  products  of  energetic  double 
and  triple-base  propellant  fall  into  the  oxidation  regime  according  to  this 
criterion. 


A  preliminary  test  matrix  was  devised  to  examine  oxidation  effects 
on  4340  barrel  steel,  the  primary  material  used  in  the  manufacture  of  gun 
barrels.  The  combination  of  test  configuration  and  gas  conditions  required 
to  achieve  the  melting  threshold  for  4340  steel  was  to  be  established  during 
the  initial  test  series.  This  threshold  condition  was  not  achieved  previously 
and  a  test  configuration  with  increased  heating, namely  the  aforementioned 
reduction  in  flow  passage  diameter,  was  instituted.  A  requirement  for  these 
and  subsequent  tests  was  obtaining  gas  pressures  and  temperatures  of  the 
magnitude  generally  obtained  in  large  caliber  weapons.  A  fundamental  property 
of  the  test  gas  mixture  which  governs  the  pressure/temperatures  produced 
in  the  shock  tube  compression  process  is  its  ratio  of  specific  heats,  Ym 
The  previous  study  concluded  that  the  optimum  Ym  for  the  facility  was  1.53. 
Therefore  the  initial  "inert"  tests  were  conducted  with  a  45.5%  nitrogen 


3.  F.A.  Vassal lo,  "Mathematical  Models  and  Congmter  Routines  Used  in 
Evaluation  of  Caseless  Aanunition  Heat  Transfer,"  Calspan  Report  No. 
GM-2948-Z-1,  June  1971. 

4.  ASM  Committee  on  Furnace  Atmospheres,  "Furnace  Atmosphere  and  Carbon 
Control,"  ASM  1964,  pg.  2. 
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54.4%  argon  mix  yielding  a  1.5  specific  heat  ratio,  and  the  reactive  gas 
mixture  tests  that  followed  were  designed  to  be  within  similar  proximity 
of  the  optimum  value. 

The  second  test  series  was  intended  to  quantify  the  contribution 
of  oxidation  by  varying  the  oxygen  concentration  at  the  4340  melting  threshold 
gas  conditions  and  by  varying  gas  conditions  with  a  constant  oxygen  concen¬ 
tration.  These  tests  were  designed  to  determine  erosion  augmentation  as 
a  function  of  oxygen  concentration,  determine  whether  or  not  melting  is  a 
prior  condition  for  increased  erosivity,  and  determine  the  effective  heat 
addition  from  oxidation.  The  basic  philosophy  behind  these  tests  was  to 
make  only  small  changes  in  a  basic  gas  mixture.  In  this  way,  test  conditions 
would  be  nearly  constant  so  that  the  true  oxidation  influence  would  be 
readily  apparent.  Furthermore,  small  increments  in  oxygen  made  it  possible 
to  generate  the  empirical  results  that  are  useful  in  guns  while  large  concen¬ 
trations  create  spectacular  but  less  useful  results. 

A  third  series  replaced  oxygen  with  carbon  dioxide.  These  tests 
were  designed  to  explore  oxidation  as  it  actually  occurs  in  guns.  Dissociated 
carbon  dioxide  is  believed  to  provide  a  major  source  of  oxygen  for  reaction 
with  barrel  steel  during  the  ballistic  cycle.  Therefore,  carbon  dioxide  was 
introduced  into  the  test  gas  in  the  amount  normally  present  in  propellant 
gas.  Gas  temperatures  were  varied  over  the  range  experienced  by  various 
propellants  of  interest.  The  object  of  these  experiments  was  to  determine 
whether  or  not  an  atmosphere  rich  in  carbon  dioxide  would  enhance  barrel 
erosion  above  the  level  attributed  to  pure  melting,  and  to  quantify  the 
erosion  in  terms  of  an  equivalent  amount  of  oxygen  for  a  given  concentration 
of  carbon  dioxide,  test  conditions  and  test  material.  This  is  the  first  step 
in  relating  chemical  erosion  in  a  gun  to  laboratory  conditions. 

Future  research  could  make  use  of  gas  mixtures  containing  both 
carbon  dioxide  and  carbon  monoxide  to  more  closely  investigate  the  neutral 
point  region,  i.e.,  C0/C0_  ~3,  where  oxidation  may  give  way  to  other  erosive 
thermochemical  reactions.  A  study  of  CO/CO.  ratios  in  excess  of  the  neutral 
point  logically  extends  into  a  subsequent  program  dealing  with  carburization. 


Table  2  .  Physical  Properties  of  4340  Steel 


Density  (p) 

Specific  Heat  (c) 

Thermal  Conductivity  (k) 
Solidus/ Liquidus 
Latent  Heat 


7850  kg/m3 
181  J/kg-°K 
10.47  J/m-sec-°K 
1720/1770  *  K 
232,600  J/kg 


(490  lbm/ft3) 

(0.14  Btu/lbm-°R) 
(19.6  Btu/ft-hr-°R) 
(2640/2720  °F) 

(100  Btu/lbm) 
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IV. 


TEST  PROCEDURE  AND  RESULTS 


4 . 1  Test  Procedures 

Collection  of  test  data  regarding  ballistics,  heating,  and  erosion 
followed  a  set  procedure.  First,  all  important  components  of  the  facility 
were  inspected  for  attrition  due  to  the  previous  firing.  Expendable  items 
such  as  "0"  rings  and  other  seals  were  replaced.  The  driven  tube,  chamber, 
and  piston  were  carefully  cleaned  to  avoid  the  presence  of  contaminants. 

The  piston  was  then  fixed  in  position  at  the  upstream  end  of  the  driven  tube. 
The  projectile  was  inserted  into  the  barrel. 

Specimens  were  characterized  prior  to  testing  using  an  analytical 
balance  for  initial  mass,  and  then  installed  within  the  chamber.  M-ll  pres¬ 
sure  gauges  were  also  installed  in  the  chamber  for  a  redundant  measure  of 
pressure.  All  required  instrumentation  including  pressure  transducers, 
thermocouples  and  velocity  screens  were  connected  to  suitable  recording 
devices:  the  piezoelectric  pressure  output  was  recorded  on  a  Bell  and  Howell 
tape  recorder  at  a  speed  of  .76  meters  per  second  (30  inches  per  second)  to 
minimize  tape  flutter  and  pulse  response  rise  time;  thermocouple  output  was 
recorded  directly  on  a  CEC  oscillograph  at  a  paper  speed  of  5  inches  per 
second;  the  projectile  transit  time  between  velocity  screens  was  recorded  by 
use  of  a  Tektronix  oscilloscope  for  1  millisecond  per  division.  The  Tektronix 
scope  was  also  used  to  record  the  pressure  signal  playback  from  the  tape  which 
was  performed  at  a  tape  speed  of  .76  meters  per  second,  the  same  speed  used 
for  recording.  Matching  the  playback  and  recording  speeds  eliminated  hys- 
teretic  pulsing  of  the  tape  heads. 

After  installation  of  the  projectile  and  specimens,  the  entire 
driver  tube  including  the  chamber  was  evacuated  to  a  pressure  level  of  less 
than  2.0  millimeters  of  mercury.  The  chamber  was  then  purged  by  filling  with 
pure  argon.  It  was  again  evacuated  and  filled  to  the  required  partial  pres¬ 
sures  with  the  gases  selected  for  the  test  mix.  These  partial  pressures  are 
dependent  upon  the  mix  ratio  desired. 

Equations  for  establishing  partial  pressure  settings  were  derived 
from  the  Dalton  model  of  ideal  gas  mixtures,  which  assumes  the  following: 

1.  The  moles  of  mixture,  n,  equals  the  sum  of  the  moles 
of  the  component  gases,  nA  ♦  nN  ♦  ng,  where  A,  N  and 
0  are  subscripts  referring  to  argon,  nitrogen  and 
oxygen,  respectively. 

2.  Each  component  gas  in  the  mixture  occupies  the  entire 
mixture  volume,  V,  which  in  this  case  is  the  volume 
of  the  driven  tube. 
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3. 


The  temperature ,  T,  of  the  components  before 
and  after  mixing  remains  constant. 


4.  The  mixture  pressure,  P,  in  this  case,  1  atmos¬ 
phere,  is  reasonably  low,  to  assure  near  ideal 


gas  behavior. 

the  components: 

PAV 

*  nART 

pNv 

*  n^RT 

pov 

*  no*T 

the  mixture: 

PV  = 

‘  nRT 

Since  V/RT  is  a  constant  in  all  the  equations: 


P 

n 


A 

A 


P 

n 


Rewriting: 


P  n 


That  is,  for  each  component  of  a  mixture  of  ideal  gases,  the  mole  fraction  and 
the  ratio  of  the  partial  pressure  to  the  total  pressure  are  equal. 

Following  these  test  preparations,  the  driver  tube  was  pressurized 
to  the  desired  level,  recording  devices  were  activated  and  the  piston  was 
released. 


After  exhausting  residual  driver  pressure,  specimens  were  removed 
and  hard  copy  was  made  of  the  test  data.  Specimens  were  inspected,  weighed 
and  measured  diametrally  at  the  aforementioned  axial  locations.  Photographs 
were  taken  of  a  sample  when  deemed  appropriate. 
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4.2 


Text  Matrix 


A  total  of  31  tests  were  conducted  during  this  program.  A  summary 
of  the  test  conditions,  the  weight  loss,  the  recession  measurements  along  the 
length  of  the  test  sample  are  given  in  Table  3.  The  initial  driver  pressure 
required  to  achieve  the  test  conditions,  the  peak  pressure,  projectile  velocity 
that  was  obtained  during  the  test  and  bulk  heating,  as  measured  by  the  in-wall 
thermocouples,  are  presented  in  Table  4. 

The  test  matrix  was  divided  into  essentially  eight  areas.  The 
first  area  involves  establishing  a  baseline  with  inert  tests.  The  test  gas 
composition  for  these  tests  was  45  1/2  percent  nitrogen  and  54  1/2  percent 
argon.  This  mixture,  when  compressed  to  a  pressure  of  275  MPa,  generated  a 
test  gas  temperature  in  the  vicinity  of  3700®K.  This  gas  temperature  is  quite 
similar  to  that  of  a  hot  double-base  propellant  such  as  M8,  After  establish¬ 
ment  of  the  inert  baseline  conditions,  three  series  of  tests  were  conducted 
with  small  concentrations  of  oxygen  to  define  the  effects  of  oxidation  in 
erosion  of  barrel  steel.  During  these  tests,  the  concentrations  of  oxygen 
were  maintained  at  a  low  level  to  minimize  the  impact  of  gas  mixture  variation 
on  gas  temperature  and  pressure.  With  this  approach,  the  effects  that  were 
observed  during  these  tests,  when  compared  to  the  inert  baseline,  were  due 
to  solely  the  oxygen  content  rather  than  changes  in  test  conditions.  This 
permitted  direct  assessment  of  the  impact  of  oxidation.  Gas  mixtures  with 
oxygen  levels  of  1,  2.5  and  5  percent  of  the  total  mixture  were  incorporated 
during  this  phase  of  the  test  program. 

Following  oxygen  tests,  carbon  dioxide  was  used  as  a  gas  consti¬ 
tuent  in  place  of  oxygen,  once  at  5  percent  concentration  and  the  remaining 
nine  runs  at  10  percent.  Variations  of  the  fraction  of  argon  and  nitrogen 
were  made  to  achieve  a  variety  of  test  gas  temperatures,  pressures,  and 
equilibrium  gas  constituents.  The  dissociation  of  carbon  dioxide  creates 
free  oxygen  at  the  higher  gas  temperatures,  which  can  promote  oxidation-type 
reactions.  The  amount  of  oxygen  available  at  peak  conditions  from  dissocia¬ 
tion  of  carbon  dioxide  is  used  as  a  means  to  interpret  the  erosion  due  to 
carbon  dioxide  in  terms  of  the  oxygen  content  of  the  previous  runs  where 
oxygen  was  introduced  as  a  primary  constituent. 

Finally,  two  short  series  of  tests  were  conducted  whereby  concen¬ 
trations  of  5  and  10  percent  hydrogen  were  introduced  with  10  percent  carbon 
dioxide.  This  created  a  gas  condition  more  nearly  like  that  in  propellant 
gases,  and  provided  metal lographic  data  for  comparison  of  the  "cracking" 
potential  of  CO2TH2  and  straight  CO2  gas  mixtures. 

Adiabatic  compression  of  the  test  gas  is  the  mechanism  by  which 
the  temperature  and  pressure  profile  is  generated.  The  test  gas  composition 
alters  the  relationship  between  the  gas  temperature  and  pressure  for  a  given 
compression  ratio.  This  relationship  is  shown  in  Figure  5  for  all  of  the  gas 
compositions  used  during  this  series  of  tests .  The  individual  points  for  the 
runs  are  shown  on  Figure  6,  as  well  as  the  previous  figure,  so  that  the  level 
of  temperature  and  pressure  is  indicated  as  well  as  the  potential  latitude 
that  can  be  achieved  in  the  Shock  Tube  Gun. 


TABLE  3.  OXIDATION  TEST  RESULTS  --  SPECIMEN  CHARACTERISTICS 
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Parenthesis  indicates  weight  gain,  or  material  buildup. 
'The  remaining  percentage  of  the  mixture  is  argon. 


TABLE  4.  HEATING/BALLISTIC  TEST  RESULTS  -  SHOCK  TUBE  GUN 
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Dashed  lines  indicate  lost  data  due  to  instrumentation  malfunction. 


Test  conditions  along  the  isentrope  were  computed  by  an  equilibrium 
combustion  program.  Appendix  II  contains  a  brief  description  of  the  code  and 
a  sample  printout.  This  program  was  also  used  to  determine  the  equilibrium 
composition  of  the  gas  mixture.  In  this  way,  with  a  specified  gas  mixture 
and  measured  pressure,  the  temperature  and  gas  constituents  generated  as  a 
result  of  the  dissociation  and  chemical  reactions  were  determined  for 
the  compression  cycle. 

The  equilibrium -gas -mixture  mole  fractions  of  the  most  significant 
constituents  are  shown  in  Table  5.  Table  5  lists  the  active  constituents 
by  mole  fraction  of  each  test  gas  mixture,  both  at  ambient,  pre-test  conditions 
and  at  the  "equilibrium,"  maximum  pressure  and  temperature  test  conditions. 
Nitrogen,  though  classified  as  "inert"  when  not  in  the  presence  of  oxygen  or 
carbon  dioxide,  has  been  included  in  the  Table  to  illustrate  its  potential 
for  forming  oxides,  acids,  and  cyanides  with  active  mixture  components 
during  a  test.  These  nitrogen  compounds,  all  gaseous  byproducts,  contribute 
little  to  direct  barrel  wear  and  erosion.  The  significant  role  they  do  play 
is  in  limiting  the  oxidizing  potential  of  the  truly  reactive  components  in 
the  mixture.  Argon  has  been  excluded  from  the  table  since  it  truly  is  "inert." 
However,  since  the  mole  fraction  of  Argon  was  not  constant  for  all  tests, 
the  mole  fractions  of  the  remaining  constituents  have  not  been  normalized  to 
equal  one.  This  practice  would  be  misleading,  particularly  in  the  CO2  tests 
where  the  total  mole  fraction  of  active  constituents  was,  for  the  majority  of 
runs,  lower  than  in  the  oxygen  tests.  In  particular,  the  mole  fraction  of 
oxygen  molecules  is  given  on  this  table.  It  is  significant  that  during  tests 
with  CO2  and  nitrogen  as  the  main  reacting  gas  species,  significant  oxygen 
was  developed.  However,  with  the  addition  of  hydrogen  into  the  mixture,  oxygen 
was  consumed  by  formation  of  water  and  the  amount  of  free  oxygen  available  in 
the  stream  at  peak  conditions  was  reduced  substantially.  This  has  a  significant 
impact  on  the  overall  results  of  this  program  which  will  be  discussed  in 
Section  V. 
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V.  DISCUSSION  AND  CORRELATION  OF  RESULTS 

5.1  Correlation  of  Erosion  Data 

5.1.1  Oxidation  Effects 

A  series  of  tests  was  conducted  to  quantify  erosion  due  to  oxi¬ 
dation  of  4340  steel  as  described  previously.  The  primary  variables  in  these 
tests  were  oxygen  concentration  of  the  test  gas  and  the  compression  ratio  which 
determined  pressure  and  temperature.  A  baseline  test  mixture  of  54  1/2  percent 
argon  and  45  1/2  percent  nitrogen  was  selected  during  the  previous  program. 1 
This  baseline  mixture  was  modified  to  include  oxygen  by  replacing  nitrogen  with 
oxygen.  The  total  oxygen  concentration  for  the  test  was  kept  relatively  low  in 
order  to  maintain  nearly  constant  gas  temperature  and  pressure  conditions  for 
the  different  concentrations  and  also  to  keep  the  net  oxidation  level  at  an 
order  of  magnitude  that  might  be  expected  in  a  gun.  The  test  series  included 
oxygen  concentrations  of  0,  1  percent,  2  1/2  percent,  and  5  percent,  and  pres¬ 
sure  levels  from  100  MPa  to  a  nominal  maximum  level  of  280  MPa. 

The  role  of  oxidation  on  surface  erosion  is  readily  shown  in 
Figure  7  where  weight  loss  is  shown  as  a  function  of  initial  oxygen  concentra¬ 
tion  at  the  maximum  nominal  pressure  of  280  MPa.  The  convective  heating  is 
sufficiently  high  to  cause  significant  material  loss  by  pure  melting  with  no 
oxygen.  Weight  loss  increases  rapidly  from  this  initial  level  with  increasing 
oxygen  concentration.  The  nearly  linear  dependence  of  weight  loss  on  oxygen 
content  indicates  a  situation  where  the  rate  of  oxidation  is  diffusion  controlled 
That  is,  the  oxidation  is  limited  by  the  number  of  oxygen  molecules  coming  in 
contact  with  the  sample  surface  rather  than  the  rate  at  which  the  reaction  occurs 

The  variation  of  weight  loss  with  flow  conditions,  as  specified  by 
pressure,  for  gas  mixtures  containing  zero  and  2  1/2  percent  oxygen  is  shown  in 
Figure  8.  Additional  inert  data  presented  in  this  figure  and  others  weTe 
obtained  during  testing  subsequent  to  this  program.5  The  threshold  at  which 
sample  weight  loss  is  first  observed  occurs  at  194  MPa  for  a  gas  mixture  con¬ 
taining  2  1/2  percent  oxygen  and  approximately  250  MPa  for  a  mixture  with  no 
oxygen.  Weight  loss  increases  as  pressure  increases  beyond  the  threshold  level. 
At  pressures  in  the  vicinity  of  and  below  the  threshold  level,  samples  tested 
with  the  2  1/2  percent  oxygen  mixture  exhibited  little  weight  change  but  the 
diametral  surface  measurements  indicated  an  oxide  buildup.  The  oxide  buildup 
at  low  pressures  gives  way  to  surface  recession  at  high  pressures. 

The  shift  in  the  erosion  threshold  is  an  indication  that  melting 
is  not  a  prior  requisite  for  the  onset  of  oxidation  and  that  surface  reactions 
are  the  mechanism  by  which  oxidation  of  steel  enhances  material  loss.  The 
erosion  versus  pressure  curves  shown  in  Figure  8  are  roughly  parallel  to  each 
other  which  indicates  that  the  effects  of  oxidation  represent  a  constant  input 
rather  than  a  driving  force  that  would  further  enhance  erosion.  The  situation 


5.  "Investigation  of  the  Role  of  Carburization  in  Gun  Barrel  Erosion  and 

Cracking,"  Bimonthly  Progress  Report  No.  3,  ARRADC0M  Contract  DAAK11-79-C 
0049,  November  1979. 
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WEIGHT  LOSS 


Figure  7  VARIATION  OF  WEIGHT  LOSS  AND  HEAT  INPUT  AS  A  FUNCTION  OF 
OXYGEN  CONCENTRATION 
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Figure  8  WEIGHT  LOSS  AND  HEAT  INPUT  FOR  VARIABLE  TEST  CONDITIONS  AND 
2-1/2%  OXYGEN  CONCENTRATION 
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where  combustion  in  the  boundary  layer  causes  catastrophic  material  erosion 
was  observed  in  the  case  of  aluminum^! 7  which  has  roughly  an  order-of-magnitude 
greater  heat  of  combustion.  In  this  instance,  catastrophic  erosion  resulted 
after  initial  melting  took  place.  Such  does  not  appear  to  be  the  case  with 
4340  steel,  but  rather,  the  surface  reaction  tends  to  serve  as  a  constant  heat 
flux  increment  to  the  convective  term. 

5.1.2  Erosion  With  Mixtures  Containing  Carbon  Dioxide 

Oxidation  in  combustion  gases  is  believed  to  occur  when  the 
ratio  of  carbon  monoxide  to  carbon  dioxide  is  less  than  the  3.0?  A  series 
of  tests  was  conducted  in  this  program  whereby  a  gas  mixture  consisting  of 
nitrogen,  carbon  dioxide  and  argon  was  compressed  to  temperatures  and 
pressures  experienced  in  guns.  The  equilibrium  concentration  ratio  of  carbon 
monoxide  to  carbon  dioxide  was  considerably  less  than  3.0  for  these  runs. 

In  many  cases,  the  equilibrium  oxygen  concentration  at  peak  pressure  was  at 
levels  similar  to  those  of  tests  where  oxygen  was  introduced  initially  as  a 
gas  constituent. 

The  initial  tests  involving  carbon  dioxide  were  conducted 
under  the  same  philosophy  as  those  involving  oxygen  with  only  a  small  concen¬ 
tration  of  carbon  dioxide  being  introduced  in  order  to  cause  only  a  small 
perturbation  to  the  test  conditions.  However,  very  little  erosion  was 
observed  during  these  initial  tests,  and  also  it  was  found  that  carbon 
dioxide,  being  significantly  different  from  nitrogen,  caused  significant 
changes  in  the  relationship  between  gas  temperature  and  pressure  as  shown 
previously  in  Figure  5.  Therefore,  the  initial  concentration  of  argon  was 
increased  in  order  to  increase  the  gas  temperature  level  and  to  increase 
the  level  of  material  erosion.  This  change  in  test  conditions  made  pressure 
alone  an  unsuitable  parameter  for  correlating  all  erosion  data. 

Convective  hot  wall  heat  flux  is  the  driving  potential  for 
melting  erosion  once  the  melting  tenperature  is  achieved.  Thus,  it  follows 
that  convective  heat  flux  to  a  hot  surface  is  a  natural  parameter  for 
correlating  simple  melting  erosion  data  resulting  from  similar  heating 
cycles.  Similarity  of  heating  cycles  is  important  because  the  heat  flux 
history  determines  the  surface  temperature.  The  flow  conditions  for  each 
test  were  computed  with  a  combination  of  two  computer  programs  and  experimen¬ 
tally  measured  pressure.  The  model  of  the  STG  cycle  (Appendix  I)  was  used 
to  create  the  proper  pressure  profile.  Also  included  in  this  particular 
code  is  a  representation  of  the  flow  through  the  test  sample  and  calcula¬ 
tion  of  the  convective  heat  flux  using  the  empirical  equation  for  turbulent 
flow  over  a  flat  plate.  The  gas  conditions  of  temperature,  density  and 


6.  Summerfield,  et.  al.,  "Erosion  of  Aluminum  by  High  Pressure  Propellant 
Gases,"  10th  JANNAF  Combustion  Meeting,  CPIA  Publication  No.  243,  August 
1973. 

7.  Fisher  and  Cytron,  "Influence  of  Combustion  on  Aluminum  Cartridge  Case 
Burnthrough,"  10th  JANNAF  Combustion  Meeting,  CPIA  Publication  No.  243, 
August  1973. 
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velocity  are  evaluated  in  this  code.  The  local  hot  wall  heat  flux  is  computed 
and  a  running  summation  of  the  heat  input  to  the  wall  is  evaluated.  The 
computed  total  heat  input  was  compared  with  the  experimentally  determined 
value  for  the  input  gas  mixture.  A  factor  was  applied  to  the  heat  flux  cal¬ 
culation  to  bring  the  total  heat  input  into  agreement  with  the  experimentally 
determined  value.  In  this  manner,  the  instantaneous  value  of  convective 
heating,  as  determined  by  the  flow  conditions  and  exclusive  of  chemical 
heating,  is  believed  to  be  reasonably  correct  for  tests  with  inert  and 
low-oxygen-concentration  gas  mixtures. 

The  gas  temperature  as  determined  by  the  STG  model  in  its  current 
state  of  development  is  only  an  approximate  calculation.  A  more  accurate 
temperature  calculation  is  provided  by  the  equilibrium  combustion  code  that 
is  described  briefly  in  Appendix  II.  This  code  computes  an  isentropic 
compression  of  gases  beginning  with  an  arbitrary  mixture.  In  addition  to 
determining  the  temperature  and  pressure,  the  concentration  of  the  various 
chemical  species  formed  during  the  equilibrium  combustion  process  are 
also  determined.  Thus,  the  mole  fraction  of  free  oxygen  formed  by  compression 
of  gas  mixtures  containing  carbon  dioxide  is  evaluated. 

An  approximate  technique  was  used  to  correct  the  heat  flux 
calculated  by  the  STG  model  to  the  more  accurate  temperature  conditions  of 
the  equilibrium  combustion  code.  The  heat  flux  to  a  surface  is  equal  to  the 
product  of  a  coefficient  and  the  temperature  difference  between  the  gas  and 
the  surface. 


q  =  h  (T0  -  Tw) . 


For  the  heat  flux  to  a  flat  plate  in  turbulent  flow,  the  coefficient  is 
functionally  proportional  to  the  gas  density,  velocity  and  viscosity. 


h  ~ 


,  .0.8  0.2 

(pu)  y 


These  gas  parameters  can  be  expressed  in  terms  of  temperature  as  follows: 


p  ~  T_1  through  the  perfect  gas  equation  of  state 

u  ~  T^2  through  the  energy  equation, 
and  p  ~  T*/2  fr0m  molecular  transport  theory. 

Thus,  the  approximate  dependence  of  h  on  the  gas  temperature  is 
h  ~  (T_1  •  T1/2)0,8  (T1/2)0,2  =  T"’4 


The  convective  hot  wall  heat  flux  computed  by  the  STG  code  is 
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qs  ~  Tgs-°‘4  (Tgs  -  Tws) 

where  T_s  and  TwS  are  the  respective  peak  values  for  gas  and  surface  tenpera- 
tures.  “Similarly,  the  convective  heat  flux  to  a  melting  steel  surface  is 


~  T  "°*4  (T 
ge  1  ge 


-  T  ) 
wm' 


where  T  is  the  gas  temperature  calculated  by  the  equilibrium  combustion  code 
and  Twg,  is  the  solidus  temperature  of  barrel  steel,  1720°K.  The  corrected 
value  for  heat  flux,  qhw»  is 


-  T 

wm 

-  T 

ws 


This  represents  an  approximate  value  of  convective  heat  flux  to  a  melting 
surface  and  was  used  in  subsequent  correlation  of  erosion  data. 

A  graph  showing  the  correlation  of  weight  loss  with  convective 
hot  wall  heating  is  shown  in  Figure  9.  All  test  gases  containing  oxygen 
initially  are  shown  with  clear  symbols  while  gas  mixtures  containing  carbon 
dioxide  are  shown  with  filled  symbols.  The  inert  as  well  as  the  21/2  per¬ 
cent  oxygen  tests  are  fairly  well  established  curves  that  show  the  shift  in 
erosion  threshold  due  to  the  addition  of  oxygen.  It  is  noted  that  test 
sequences  for  mixtures  containing  1  percent  oxygen  and  5  percent  oxygen  are 
not  sufficiently  complete  to  define  the  onset  of  erosion,  however,  the  points 
from  those  curves  fall  in  proper  sequence  with  respect  to  the  inert  and  2  1/2 
percent  oxygen  curves  so  that  at  any  point  there  is  nearly  a  linear  increase 
in  erosion  with  increasing  oxygen  as  noted  earlier  in  Figure  7. 

The  data  for  CO2  runs  show  that,  in  some  cases,  carbon  dioxide  tends 
to  enhance  erosion,  while  in  others,  the  erosion  is  about  the  same  as  for 
the  inert  gas  mixture.  It  is  interesting  to  note  that  the  amount  of  erosion 
above  that  for  the  inert  gas  can  be  correlated  reasonably  well  with  the  amount 
of  free  oxygen  at  peak  conditions.  In  particular,  tests  50,  51  and  53  all 
contain  approximately  1  percent  or  more  free  oxygen  at  peak  conditions.  The 
points  from  these  tests  all  fall  above  the  curve  for  2  1/2  percent  oxygen. 
Tests  52  and  54,  contain  0.6  and  0.7  percent  free  oxygen,  respectively,  and 
these  points  fall  above  the  curve  for  1  percent  oxygen.  Finally,  the  weight 
loss  experienced  during  run  55,  which  contained  0.5  percent  free  oxygen  at 
peak  conditions  is  slightly  below  the  values  for  1  percent  oxygen.  Thus, 
there  appears  to  be  twice  the  erosion  for  these  runs  than  for  runs  with  an 
equivalent  amount  of  oxygen.  This  trend  was  broken  by  run  47  from  which  there 
was  no  measurable  erosion  and,  yet,  it  contained  0.3  percent  free  oxygen. 

The  point  fell  slightly  below  the  inert  curve,  very  close  to  the  threshold 
of  erosion. 
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CONVECTIVE  HOT  WALL  HEAT  FLUX  ~  J/mmz-S«c 

INFLUENCE  OF  CONVECTIVE  HEATING  AND  OXIDATION-TYPE 
REACTIONS  ON  EROSION  OF  4340  STEEL 


If  one  examines  this  figure,  it  is  apparent  that  most  of  the  tests 
with  CO2  that  experienced  erosion  in  excess  of  the  inert  case  fell  either  on 
or  very  close  to  the  line  of  constant  heat  flux  at  the  280  joules  per  square 
millimeter  per  second.  The  weight  loss  versus  equilibrium  oxygen  concentration 
at  this  heating  level  for  both  the  tests  with  oxygen  and  with  carbon  dioxide  is 
plotted  in  Figure  10.  In  this  figure,  the  augmentation  of  erosion  due  to  car¬ 
bon  dioxide  over  that  of  an  equivalent  amount  of  oxygen  becomes  readily  apparent. 
It  is  evident  from  these  two  figures  that  the  equilibrium  concentration  of 
oxygen  in  the  stream  is  an  indicator  of  erosion  enhancement  over  the  inert  case. 

The  enhanced  erosion  observed  during  tests  with  carbon  dioxide  over 
those  during  tests  with  oxygen-nitrogen  mixtures  suggests  the  existence  of 
an  erosion  mechanism  other  than  oxidation.  One  might  postulate  that  oxidation 
occurring  initially  at  the  surface  might  cause  depleted  oxygen  in  that  region 
which  would  promote  formation  of  carbides.  Carbon-rich  steels  and  carbides 
have  lower  melting  temperatures  than  pure  steels  or  oxides  as  shown  in  Figure 
11,  which  suggests  one  mechanism  by  which  enhanced  erosion  in  the  presence  of 
CO2  might  occur.  4340  steel  melts  at  approximately  1700°K,  while  iron  carbide 
is  believed  to  melt  in  the  vicinity  of  1440®K,  The  term  "melting"  is  used 
loosely  here  to  mean  the  solidus  point  or  onset  of  the  mushy  state  indicated 
in  Figure  11.  It  is  assumed  that  shear  and  deformation  by  high  pressure  gas 
flow  removes  the  material  at  this  temperature.  For  the  case  of  run  51,  just 
this  change  in  melting  temperature  could  account  for  a  13  percent  increase  in 
material  loss  rate.  As  the  carbon  content  increases  in  steel,  its  thermal 
conductivity  decreases,  so  that  for  a  given  heat  flux  to  the  surface,  one  might 
expect  accelerated  removal.  Such  mechanisms  represent  only  conjecture  at  this 
point.  Future  analyses  of  test  samples,  however,  may  indicate  an  elevated 
carbon  content  at  the  surface  which  would  support  the  argument. 

Part  of  the  enhanced  CO7  erosion  indicated  by  Figure  10  may  be  attri¬ 
buted  to  differences  in  the  ballistic  cycle.  The  amount  of  material  removal 
is  not  only  a  function  of  the  convective  hot  wall  heat  flux,  but  on  its  entire 
heating  history.  Both  the  time  at  the  high  heating  rate  and  the  dynamic 
heat  conduction  processes  are  important.  The  measured  total  heat  input 
represents  the  integrated  convective  and  chemical  heat  flux  to  the  sample 
and  is  an  indicator  of  differences  in  ballistic  cycle  time  between  tests  with 
equal  levels  of  hot  wall  heat  flux  in  a  similar  chemical  environment. 

The  measured  total  heat  input  is  given  in  Table  4.  This  represents 
the  residual  heat  in  the  sample  at  the  end  of  the  test.  This  quantity  is  in 
error  for  samples  experiencing  large  erosion  because  the  act  of  material  removal 
also  removes  heat  from  the  location  of  the  heat  sensor.  Therefore,  a  correction 
factor  was  added  to  the  heat  input  data  for  the  purpose  of  correlating  weight 
loss  measurements.  The  factor  is  based  on  the  assumptions  that  the  material 
removed  was  at  1720°K  and  that  the  material  was  sloughed  away  before  the  latent 
heat  was  absorbed.  The  correction  factor  is 

AQ  =  pcAx  AT  =  3.30  Ax  J/mm^ 

which  is  simply  added  to  the  heat  input  given  in  Table  4  to  form  Qcorr-  The 
value  of  Ax  is  the  diametral  recession  under  column  B,  which  represents  the 
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location  of  the  heat  sensor.  If  the  latent  heat  term  were  included,  the 
correction  factor,  AQ,  would  be  increased  by  0.91  Ax. 

The  correlation  of  weight  loss  with  the  corrected  heat  input  is  shown 
in  Figure  12.  With  the  exception  of  a  few  "flyers,"  the  dotted  line  connects 
the  points  in  a  reasonable  fashion.  The  weight  loss  remains  at  zero  with 
increasing  heat  input  until  the  combination  of  heat  input  and  peak  hot  wall 
heating  is  sufficiently  high  to  cause  material  removal.  It  is  noted  that 
several  samples  experienced  total  heating  in  excess  of  the  apparent  threshold 
but  did  not  lose  material.  During  these  tests,  the  ballistic  cycle  was  long 
which  gave  a  high  total  heat  input,  but  the  peak  hot  wall  heat  flux  was  below 
the  level  required  to  achieve  melting  conditions.  Of  greater  significance, 
however,  is  the  observation  that  data  from  runs  50  and  51  appear  to  have  higher 
weight  loss  than  those  tests  with  a  gas  containing  2-1/2  percent  oxygen  because 
the  total  heat  input  is  higher.  The  cause  of  the  elevated  heat  input  is  not 
clear  at  present.  Measured  pressure  curves  indicated  the  ballistic  cycle  for 
tests  with  carbon  dioxide  were  nearly  10  percent  shorter  in  duration  than  for 
the  tests  with  oxygen.  The  increased  total  heating  is  likely  a  combined  effect 
of  surface  chemistry,  convective  heat  transfer  coefficient,  and  accumulated  low 
level  heating  before  a  measurable  pressure  was  reached.  Clearly,  this  is  an 
area  that  requires  additional  investigation. 

It  has  been  shown,  however,  that  some  gas  mixtures  containing  carbon 
dioxide  are  more  erosive  than  an  inert  gas  at  similar  heating  environments  and 
that  the  level  of  erosivity  appears  to  correlate  with  the  equilibrium  level  of 
free  oxygen  in  the  gas  at  peak  conditions.  This  is  in  contrast  to  the  work  of 
Caveny8  where  no  enhanced  erosion  in  an  atmosphere  of  carbon  dioxide  without 
water  was  observed.  The  current  result  may  be  attributed  to  the  factor-of-four 
longer  test  times  of  the  STG  over  Caveny's  ballistic  compressor. 

If  hydrogen  is  added  to  the  gas ,  the  amount  of  free  oxygen  is  decreased 
substantially  due  to  formation  of  water.  In  the  situations  tested,  no  erosion 
was  observed  in  excess  of  the  inert  erosion  for  gases  containing  mixtures  of 
carbon  dioxide  and  hydrogen.  Tests  were  only  conducted  up  to  the  threshold  of 
erosion  of  an  inert  gas  mixture.  Additional  tests  are  required  to  extend  the 
data  base  to  higher  heating  conditions  for  hydrogen-carbon  dioxide  mixtures. 

5.2  Erosion  Distribution 

The  distribution  of  surface  erosion  along  the  length  of  the  test  sample 
is  shown  in  Figure  13  for  gas  mixtures  with  0,  2-1/2  percent,  and  5  percent  oxy¬ 
gen.  The  pressures  and  temperatures  were  approximately  the  same  for  each  of  the 
tests  so  differences  between  the  three  curves  are  primarily  due  to  the  oxygen 
concentration.  It  is  noted  that  the  data  for  0  percent  oxygen  show  moderately 
high  erosion  at  the  first  point  of  measurement.  The  amount  of  erosion  decreases 
rapidly  until  there  is  essentially  no  erosion  at  the  third  and  fourth  points  of 
measurement,  which  are  representative  of  the  downstream  half  of  the  test  sample. 


8.  Caveny,  L.,  "Erosion  of  Steel  by  Combustion  Gases,"  presented  at  ARO-sponsored 
Workshop  II  on  Mechanisms  of  Erosion  in  Hot  Flowing  Media,  9-12  October  1979. 
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%  EQUILIBRIUM  02  CONCENTRATION 


Figure  10  CORRELATION  OF  WEIGHT  LOSS  OF  4340  STEEL  WITH  TEST  GAS 

EQUILIBRIUM  CONCENTRATION  OF  OXYGEN  AT  A  PEAK  CONVECTIVE 
HOT  WALL  HEATING  RATE  OF  280  J/mm2-rec 
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Figure  13  INFLUENCE  OF  OXIDATION  ON  SURFACE  RECESSION 
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The  data  from  the  sample  tested  at  a  2  1/2  percent  concentration 
of  oxygen  show  a  somewhat  higher  erosion  at  the  first  measurement  point  but 
the  same  basic  distribution  over  the  length  of  the  sample  as  for  the  inert 
gas  mixture.  The  erosion  at  the  flow  exit  point  was  very  small.  For  a  5  per¬ 
cent  oxygen  mixture,  the  sample  experienced  a  heavy  erosion  over  its  entire 
length  with  the  erosion  at  the  first  three  stations  being  nearly  the  same  and 
dropping  off  only  slightly  at  the  flow  exit  point. 

These  observations  are  interpreted  with  the  help  of  Figure  14. 

Here,  characteristic  heat  flux  distributions  are  shown  plotted  versus  distance 
from  the  sample  stagnation  point.  The  distribution  of  heating  over  the  curved 
portion  of  the  sample  to  the  sonic  point,  which  is  denoted  by  S*,  is  approximated 
on  the  basis  of  melting  patterns  observed  near  the  melting  threshold.  From 
that  point,  the  sample  is  assumed  to  be  subjected  to  a  distribution  of  heat 
flux  that  is  analogous  to  the  heat  flux  to  a  flat  plate  in  turbulent  boundary 
layer  flow.  The  inert  heating  profile  was  arbitrarily  established  for 
melting  to  occur  at  the  leading  edge. 

Several  assumptions  can  be  made  as  to  the  impact  of  oxygen  on 
the  distribution  of  erosion.  It  was  observed  previously,  that  the  presence 
of  oxygen  in  the  gas  mixture  caused  a  shift  in  the  threshold  of  erosion.  This 
indicates  that  oxidation  increases  the  effective  level  of  total  heating  which 
causes  the  melting  threshold  to  occur  at  a  lower  level  of  convective  heating 
and  less  severe  flow  conditions.  At  downstream  locations  on  the  test  sample, 
the  surface-chemistry-effect  might  be  expected  to  elevate  the  heat  flux  level 
by  an  increment  as  illustrated  conceptually  in  Figure  14. 

Once  the  level  of  heating  required  to  achieve  melting  is  reached, 
the  rate  at  which  material  removal  occurs  at  a  given  location  is  greatly 
accelerated  by  any  additional  heating.  Therefore,  most  of  the  difference  in 
surface  recession  between  the  5  and  2  1/2  percent  oxygen  tests  can  be 
attributed  to  an  increase  in  chemical  heating  at  downstream  locations  from 
threshold  to  levels  substantially  in  excess  of  the  threshold.  The  percentage 
increase  in  heating  in  excess  of  that  required  for  melting  is  much  less  at 
the  forward  stations  so  changes  there  are  not  as  dramatic.. 

Melting  in  the  forward  portion  of  the  sample  and  subsequent 
combustion  results  in  addition  of  heat  to  the  boundary  layer,  which  increases 
the  effective  gas  temperature.  This  augments  the  convective  heating,  in 
addition  to  local  heating  due  to  surface  chemistry  at  downstream  locations. 

The  magnitude  of  this  effect  is  governed  by  the  heat  released  during  the 
oxidation  process  and  the  distribution  of  this  energy  in  the  stream.  This 
was  found  to  be  a  dominant  factor  for  aluminum  cartridge  case  burn-through6 
because  aluminum  has  a  high  heat  of  reaction  with  oxygen.  FeO  has  a  heat 
of  formation  only  15  percent  that  of  aluminum  oxide  and  effects  of  combus¬ 
tion  in  the  boundary  layer  are  believed  to  be  correspondingly  small. 
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Figure  14  INFLUENCE  OF  OXIDATION  ON  SURFACE  HEATING 


5.3 


Surface  Features 


The  Scanning  Electron  Microscope  (SEM)  was  used  to  examine  the 
surface  features  of  the  samples  tested  during  this  program.  Melt  layers, 
oxide  coatings,  and  cracking  are  among  those  features  that  were  observed 
during  this  examination.  The  samples  were  sliced  end  to  end  to  permit  direct 
viewing  of  the  surface.  All  samples  were  examined  in  this  manner  but  several 
that  are  representative  of  the  overall  observations  are  discussed  here. 

Sample  1  from  Run  32  was  tested  at  a  nominal  pressure  of  280  MPa 
in  a  nitrogen- argon  environment.  Internal  diameter  measurements  indicate  the 
sample  experienced  virtually  no  surface  recession  except  near  the  leading  edge 
SEM  photographs  of  the  flow  entrance,  center,  and  flow  exit  regions  of  the 
sample  are  shown  in  Figure  15 .  The  surface  near  the  flow  entrance  has  a 
grainy  texture  and  an  extremely  fine  check  pattern.  The  surface  roughness  is 
probably  related  to  original  machining  marks  and  grain  boundaries  in  the  steel 
The  center  of  the  sample  appears  to  exhibit  a  slightly  thicker  melt  layer 
speckled  with  solidified  droplets  that  originated  farther  upstream.  The  fine 
crack  structure,  which  does  not  resemble  cracking  in  a  gun,  also  persists  at 
this  location.  Finally,  it  is  evident  that  no  surface  melting  took  place  at 
the  sample  flow  exit,  as  indicated  by  the  presence  of  machining  grooves 
beneath  solidified  spray. 

Sample  Nos.  2  and  8  were  tested  at  nominal  pressures  of  280  MPa 
in  an  atmosphere  containing  1.2  percent  oxygen  and  5  percent  oxygen,  respec¬ 
tively.  These  samples  illustrate  surface  features  that  are  characteristic 
of  all  tests  conducted  at  high  pressure  in  various  concentrations  of  oxygen. 
The  SEM  photographs  of  Figures  16  and  17  at  the  entrance,  center,  and  exit 
are  similar  inasmuch  as  particles  protrude  from  what  appears  to  be  a  melt 
layer.  At  lower  magnification,  this  surface  has  a  crust-like  appearance  as 
shown  in  Figure  18 .  The  melt  layer  contains  a  fine  microcrack  pattern  of 
much  the  same  character  as  the  test  without  oxygen.  This  also  bears  little 
resemblance  to  the  cracking  that  has  been  observed  in  gun  tubes. 

A  series  of  tests  were  conducted  in  an  atmosphere  containing  a 
2  1/2  percent  concentration  of  oxygen  for  a  range  of  peak  pressures.  The 
surface  features  illustrated  in  Figures  19 ,  20  and  21  show  the  progression  of 
surface  cracking  from  test  conditions  that  represent  the  onset  of  melting 
erosion  to  those  at  the  higher  melting  erosion  conditions.  The  surface 
features  at  the  sample  inlet  appear  to  be  quite  rough  at  the  onset  of  erosion. 
This  rough  surface  feature  then  gives  way  to  what  appears  to  be  a  molten 
sublayer  with  particles  protruding  at  the  higher  pressures.  The  most 
significant  aspect  of  this  series  of  three  figures  is  the  crack  patterns 
that  are  present  at  the  center  of  the  sample.  The  cracks  that  occurred  during 
Run  34,  in  particular, are  somewhat  smaller  in  magnitude,  but  have  the  same 
general  characteristics  as  those  observed  in  guns.  The  flow  exit  region  is 
covered  with  resolidified  deposits  for  the  low  pressure  test  but  this  gives 
way  to  the  general  eroded  melted  layer  with  some  resolidifed  deposits  at  the 
higher  pressures. 


SAMPLE  1.  RUN  32,  1000  MAGNIFICATION 

Figure  15  SURFACE  FEATURES  OF  A  SAMPLE  MADE  FROM  4340  STEEL  AND  TESTED  IN 
AN  INERT  ATMOSPHERE  AT  275  MPa 


SAMPLE  2,  RUN  33,  1000  MAGNIFICATION 

Figure  16  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  5  PERCENT  OXYGEN 
ATMOSPHERE  TO  283  MPa,  3516°K 
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SAMPLE  8.  RUN  42,  1000  MAGNIFICATION 

Figure  17  SURFACE  FEATURES  OF  A  SAMPLE  MADE  FROM  4340  STEEL  AND  TESTED  IN 
AN  ATMOSPHERE  CONTAINING  1.2  PERCENT  OXYGEN  AT  275  MPa 


SAMPLE  4,  RUN  34.  1000  MAGNIFICATION 

Figure  19  SURFACE  FEATURES  OF  A  4340  STEEL  SAMPLE  TESTED  IN  AN  ATMOSPHERE 
CONTAINING  2-1/2  PERCENT  OXYGEN  AT  175  MPa 
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SAMPLE  7,  RUN  41, 1000 

Figure  20  4340  STEEL  SAMPLE  SURFACE  1 
ATMOSPHERE  TO  200  MPa,  3237' 


It  is  interesting  to  note  that  the  sample  tested  at  low  pressure 
experienced  no  weight  loss  and  a  surface  buildup.  All  samples  tested  in 
oxygen  at  pressures  from  170  to  230  MPa  experienced  surface  buildup  and 
cracking  at  the  center.  At  higher  pressures,  this  cracking  feature  seemed  to 
diminish  until  a  very  fine  crack  structure  was  observed  at  pressures  in  the 
vicinity  of  280  MPa.  It  is  postulated  that  cracking  may  be  associated  with 
oxide  thickness  and  perhaps  melting  of  the  oxide  layer.  At  the  higher  pres¬ 
sures,  melting  and  material  removal  occurred  over  the  length  of  the  sample 
but  at  the  lower  pressures  the  buildup  was  bbserved.  This  buildup,  which  is 
believed  to  be  an  oxide  layer,  appeared  to  crack  upon  cooling.  The  presence 
of  an  oxide  or  other  form  of  brittle  surface  layer  may,  indeed,  be  related  to 
initiation  of  cracking  in  gun  barrels. 

The  surface  features  of  samples  tested  in  a  10  percent  carbon 
dioxide  gas  mixture  ?  ;d  experiencing,  in  one  case,  a  high  degree  of  erosion 
and,  in  the  other  case,  virtually  no  erosion,  are  shown  in  Figures  22  and  23. 
respectively.  Test  sample  No.  14  was  tested  during  Run  48  in  a  gas  mixture 
that  contained  35  1/2  percent  nitrogen.  The  final  gas  temperature  was  approxi¬ 
mately  3300°K.  Sample  No.  17  was  tested  during  Run  51  in  a  mixture  that  con¬ 
tained  only  5  percent  nitrogen  so  that  its  gas  temperature  was  3800°K.  The 
peak  pressures  were  306  MPa  during  Run  48  and  285  MPa  during  Run  51. 

During  Run  51  the  sample  experienced  large  erosion  over  its  entire 
length  (Figure  22)  and  the  SEM  photographs  of  the  surface  give  the  appearance 
that  the  surface  was  scrubbed  and  eroded  quite  heavily.  The  cracks  that  are 
observed  are  very  small  microcracks.  Other  aspects  of  the  structure  appear 
to  be  protrusions  or  small  particles  that  were  engulfed  in  the  melt  layer. 

Surface  features  were  quite  different  for  the  sample  tested  during 
Run  48  (Figure  23)  at  the  center  and  exit  locations.  The  center  in  particular 
shows  an  interesting  phenomenon.  A  rather  thick  region  of  melt  covers  most  of 
the  observed  surface  area.  However,  in  the  center  of  the  picture  is  a  region 
where  the  melt  layer  is  either  missing  entirely  or  is  extremely  thin.  The 
surface  cracking  is  very  small  and  practically  insignificant  in  that  region 
not  covered  by  the  melt.  Cracking  and  surface  checking  are  readily  observed 
in  the  region  covered  by  the  thick  melt  layer.  This  observation  appears  to 
correlate  with  a  similar  observation  during  tests  with  oxygen. 

Finally,  Sample  26,  tested  during  Run  60  in  a  gas  mixture  that 
consisted  of  10  percent  carbon  dioxide  and  5  percent  hydrogen,  exhibited  rather 
severe  cracking  as  shown  in  Figure  24.  There  was  virtually  no  material  removed 
during  this  test.  A  slight  buildup  was  observed  over  a  portion  of  the  sample 
but  the  cracks  observed  here  are  quite  similar  to  those  observed  previously 
during  tests  involving  oxygen  near  the  threshold  of  erosion  when  a  surface 
buildup  was  recorded.  The  surface  of  this  sample  is  somewhat  smoother  but 
otherwise  quite  similar  to  those  tested  in  mixture  involving  oxygen  and  carbon 
dioxide  without  hydrogen. 
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SAMPLE  14,  RUN  48,  1000  MAGNIFICATION 
Figure  23  SURFACE  FEATURES  OF  A  4340  STEEL  SAMPLE  TESTED  IN  AN 
ATMOSPHERE  CONTAINING  10  PERCENT  CARBON  DIOXIDE  AND 
35-1/2  PERCENT  NITROGEN  AT  306  MPa 
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SAMPLE  26,  RUN  60,  1000  MAGNIFICATION 

Figure  24  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  10  PERCENT  CARBON  DIOXIDE 
-5  PERCENT  HYDROGEN  ATMOSPHERE  TO  210  MPa,  3316°K 


5  .-1  Metal lographic  Examination 

The  test  samples  were  sectioned,  mounted,  and  polished  so  that  the 
material  might  be  examined  in-depth  at  the  location  of  the  in-wall  thermocouple 
After  polishing,  each  of  the  test  samples  was  etched  in  a  2  percent  Nital  solu¬ 
tion.  This  is  the  procedure  normally  used  to  observe  the  white  layer  on  the 
surface  of  a  gun  barrel  and  the  thermally  altered  layer,  sometimes  referred 
to  as  the  martensite  or  hard  layer  that  extends  into  the  barrel  material. 

White  layers  in  gun  barrels  have  been  the  subject  of  extensive 
investigations  and  studies.®  It  has  been  stated^  that  white  layers  are  formed 
primarily  in  reducing  atmospheres  containing  carbon  dioxide  and  carbon  monoxide 
During  the  current  tests,  however,  white  layers  were  observed  in  gas  mix¬ 
tures  of  oxygen,  nitrogen,  and  argon  in  the  absence  of  CO  and  CO2.  The 
white  layers  were  observed  during  tests  where  there  was  a  material  buildup 
on  the  surface  after  the  test  as  well  as  those  where  material  was  removed 
during  the  test.  An  example  of  a  white  layer  formed  in  the  presence  of 
oxygen  is  shown  in  Figure  25. 

White  layers  as  distinct  as  those  cited  above  were  not  observed 
during  tests  with  carbon  dioxide.  An  example  of  the  thermally  altered  layer 
with  a  practically  non-existent  white  layer  is  shown  in  Figure  26.  In  this 
figure,  the  white  layer  is  more  prominent  than  in  any  other  photograph  made 
with  a  carbon  dioxide  or  carbon  dioxide-hydrogen  gas  mixture.  It  is  also 
noted  that  the  white  layer  was  not  observed  during  a  test  with  a  gas  mixture 
of  nitrogen  and  argon  as  illustrated  in  Figure  27.  This  seems  to  indicate 
that,  during  these  tests,  the  formation  of  oxide  may  have  been  the  source  of 
the  white  layer  as  opposed  to  a  carbon  enriched  layer.  The  formation  of  this 
layer  is  not  fully  understood  at  this  point  and  additional  research  is  required 
to  identify  other  causes  for  the  formation  of  the  white  layers  as  they  are 
observed  in  guns. 

The  thermally  altered  layer  that  extends  beneath  the  white  layer 
is  caused  by  heating  the  material  above  the  temperature  required  to  form  the 
austentitic  reaction  (  ~1150°K)  and  then  quenching  it  sufficiently  rapidly 
to  form  hard  martensite.  The  thickness  of  this  layer,  which  indicates  the 
depth  at  which  these  conditions  were  satisfied,  is  illustrated  in  Figure  28. 
Also  shown  are  calculated  temperature  distributions  into  the  sample  for 
different  times  past  the  peak  heating  condition.  In  this  curve,  it  is  shown 
that  the  depth  of  the  thermally  altered  layer  corresponds  approximately  to 
the  level  of  1150°K  and  the  point  at  which  this  level  is  maintained  for  an 
extended  period  of  time  as  a  result  of  the  heating  and  quenching  action. 
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SAMPLE  11,  RUN  45,  200 MAGNIFICATION 


Figure  25  WHITE  AND  THERMALLY  ALTERED  LAYERS  CREATED  BY  A  TEST  GAS 
CONTAINING  2-1/2  PERCENT  OXYGEN  AND  AT  A  PRESSURE  OF  254  MPa 


SAMPLE  20,  RUN  54,  400  MAGNIFICATION 

Figure  26  THERMALLY  ALTERED  LAYER  CREATED  BY  A  TEST  GAS  CONTAINING 
10  PERCENT  CARBON  DIOXIDE  AT  A  PRESSURE  OF  308  MPa.  NOTE  THE 
ABSENCE  OF  A  PROMINENT  WHITE  LAYER 


SAMPLE  1,  RUN  32,  400  MAGNIFICATION 

Figure  27  THERMALLY  ALTERED  LAYER  CREATED  BY  A  TEST  GAS  CONTAINING 

45.5  PERCENT  NITROGEN  AT  A  PRESSURE  OF  281  MPa.  NOTE  THE  ABSENCE 
OF  A  WHITE  LAYER. 


Figure  28  RELATIONSHIP  BETWEEN  THE  HARD  LAYER  AND  THE  TEMPERATURE  PROFILE 


VI.  CONCLUSIONS 


This  program  is  part  of  an  ongoing  effort  to  isolate  the  various 
phenomena  that  contribute  to  gun  barrel  erosion.  Erosion  of  4340  steel  in 
an  oxidizing  atmosphere  was  the  primary  topic  of  investigation.  This 
research  led  to  the  following  conclusions  and  accomplishments: 

1.  Addition  of  small  quantities  of  oxygen  to  the  gas  mixture 
was  shown  to  increase  the  erosion  of  4340  steel. 

2.  The  increase  in  erosion  was  nearly  linear  with  oxygen 
content  which  is  indicative  of  diffusion  controlled 
oxidation. 

3.  The  addition  of  oxygen  shifted  the  erosion  threshold  to 
less  severe  flow  conditions. 

4.  An  inert  gas  mixture  produced  no  indication  of  surface 
melting  at  conditions  comparable  to  the  erosion  threshold 
for  a  gas  mixture  containing  2  1/2  percent  oxygen.  This 
indicates  that  oxidation  takes  the  form  of  a  surface  reac¬ 
tion  rather  than  combustion  of  molten  droplets. 

5.  Erosion  in  excess  of  that  for  an  inert  gas  mixture  was 
observed  with  gas  mixtures  containing  carbon  dioxide  at 
sufficiently  severe  flow  conditions  and  where  the  CO/CO2 
ratio  was  much  less  than  3.0. 

6.  The  amount  of  erosion  in  the  CO2  atmosphere  correlated  with 
the  calculated  amount  of  oxygen  at  peak  conditions  assuming 
that  the  gas  was  in  chemical  equilibrium. 

7.  The  correlation  of  erosion  with  equilibrium  oxygen 
indicated  that,  in  most  cases,  erosion  was  more  severe 
in  the  gas  containing  carbon  dioxide  than  that  initially 
mixed  with  oxygen.  The  reason  for  the  apparent  increase 
in  heating  is  not  clear  at  present. 

8.  The  erosion  observed  with  carbon  dioxide  may  be  applicable 
to  hot  double  base  propellants  such  as  M8. 

9.  Cracks  were  observed  during  virtually  all  tests  except 
those  with  inert  gases  and  were  most  prominent  near  the 
threshold  of  erosion  where  surface  layers  were  thickest. 

10.  White  layers  were  observed  in  an  oxidizing  atmosphere  in 
contrast  to  the  findings  of  Reference  10. 
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APPENDIX  I 


SHOCK  TUBE  GUN  COMPUTER  SIMULATION 


Overview  | 

The  mathematical  model  described  in  this  Appendix  simulates  the  j 

operation  of  the  Shock  Tube  Gun.  The  simulation  provides  a  complete  descrip¬ 
tion  of  the  entire  cycle  of  the  Shock  Tube  Gun,  beginning  with  release  and  J 

subsequent  acceleration  of  the  piston  by  the  high  pressure  driver  gas.  As  j 

the  piston  is  accelerated  through  the  driven  tube,  the  simulation  computes 
the  increase  in  pressure  and  temperature  of  the  test  gas.  In  addition,  the 
simulation  evaluates  the  total  temperature  pressure  and  density  in  the 
plenum  chamber  and  computes  the  flow  through  the  test  specimen.  The  heat 
flux  to  the  specimen,  the  resulting  temperature  history  at  a  location  on 
the  surface  of  the  specimen  and  the  temperature  distribution  normal  to  I 

the  surface  are  also  calculated.  Finally,  the  simulation  calculates  the  J 

travel  of  a  projectile  through  the  barrel.  ! 

|  : 

The  objective  of  this  code  is  to  provide  a  means  for  calculating 

test  conditions  for  the  purpose  of  establishing  the  initial  driver  pressure  | 

and  gas  mixture.  Differences  in  the  ballistic  cycle  of  the  Shock  Tube  Gun  ‘  j 

due  to  gas  composition  are  reflected  through  differences  in  pressure  and  j 

total  heat  input.  A  primary  use  of  the  code  is  to  help  distinguish  between  ;  ] 

erosion  due  to  melting  and  that  due  to  chemical  effects.  This  is  done 

through  computation  of  the  convective  heating  to  the  sample  without  chemistry  j 

which  provides  a  means  for  comparing  tests  within  a  test  matrix  on  an  equal 

basis  as  far  as  the  inherent  flow  heating  of  the  test  gas.  Thus,  excess 

material  removal  from  one  gas  mixture  in  comparison  to  another  is  likely  j 

due  to  chemical  effects.  This  further  allows  estimates  to  be  made  of  the 

effective  heat  input  due  to  chemical  effects.  This  can  be  done  by  comparing  j 

the  heat  input  at  the  onset  of  erosion,  or  at  points  of  equal  erosion  between 

inert  and  chemically  active  gases.  The  code  enables  the  facility  user  to  j 

quantify  levels  of  heating  experienced.  i  \ 

;  \ 
i  i 

The  major  assumption  applied  with  formulating  the  code  was  that  I 

of  quasi-steady  operation.  That  is,  pressure  waves  and  other  unsteady  aspects 
of  the  event  are  not  calculated.  The  pressure  is  assumed  to  be  constant 
throughout  the  driver  system  and  throughout  the  driven  tube  at  any  instant 
of  time  during  the  compression  cycle.  ! 

i 

The  other  limiting  assumption  that  is  currently  employed  in  this  j 

code  is  that  of  a  frozen  gas  composition  whereby  the  initial  gas  composition  j 

is  assumed  to  be  maintained  throughout  the  ballistic  cycle.  This  assumption 
influences  the  resulting  temperature  and  pressures  to  some  extent  in  cases  where 
chemical  reactions  and  dissociation  become  important. 

The  individual  gas  constituents  are  assumed  to  be  mixed  uniformly. 
Furthermore,  the  parameters  used  by  the  Van  der  Waal  equation  of  state  for  j 


an  imperfect  gas  and  the  temperature  dependency  of  specific  heat  are  assumed 
to  be  satisfied  through  a  linear  averaging  according  to  mole  fraction. 


The  assumption  of  frozen  gas  constituency  is  an  interim  assumption 
that  will  be  relieved  when  a  chemical  equilibrium  gas  code,  also  discussed 
in  this  appendix,  can  be  combined  with  the  STG  model.  At  that  time,  the  gas 
constituency  will  be  assumed  to  be  in  thermal  and  chemical  equilibrium  at 
all  times. 


Model  Description 

Piston  Motion 

The  piston  motion  is  evaluated  by  applying  a  force  balance  on  the 
piston.  The  accelerating  force  is  applied  by  the  high-pressure  nitrogen 
driver,  500  to  700  psi,  on  the  upstream  side  of  the  piston.  The  driven  gas 
on  the  downstream  side  of  the  piston  is  initially  at  atmospheric  pressure. 

It  is  assumed  that  no  gas  leaks  past  the  piston.  This  is  essentially  verified, 
at  least  initially,  by  the  maintenance  of  a  perfect  seal  and  the  ability  to 
evacuate  the  driven  gas  chamber.  The  nitrogen  gas  is  represented  by  the 
ideal  gas  equation  of  state.  The  driver  conditions  are  calculated  from 
conservation  of  energy  principles  which  are  used  to  continually  evaluate  the 
amount  of  energy  that  is  being  transferred  from  the  gas  to  the  piston.  As 
cited  previously,  this  calculation  is  quasi-steady  in  that  the  unsteady 
expansion  aspects  are  not  considered  and  the  pressure  is  assumed  to  be 
constant  throughout  the  driver  system. 

The  driver  gas  (nitrogen)  properties  are  assumed  to  be  constant 
over  the  range  of  temperature  and  pressure  encountered  during  the  compress ion 
cycle  and  are  specified  by: 

Equation  of  state  gas  constant,  R  =  55.0  ft-lbf/lbm°R 
Specific  heat  at  constant  volume,  Cy  =  0.177  Btu/lbm-°R 
Specific  heat  at  constant  pressure,  Cp  =  0.248,  and 
Ratio  of  specific  heats,  y  =  1.4 

The  assumption  of  a  perfect  seal  at  the  piston  infers  the  existence 
of  constant  gas  mass  in  the  driver  system  during  the  cycle  so  that 


where  R  and  m  are  constants,  the  driver  volume,  V,  is  expressed  in  terms  of 
piston  travel  and  the  initial  volume  by 

v  -  v0  +  Vp, 
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where  Aj  is  the  driver  tube  area. 

Gas  temperature  is  expressed  in  terms  of  internal  energy  where 
E  =  CylllpT 

and  the  energy  change  resulting  from  work  expended  through  piston  motion 
is : 

AE  =  pdAtAxp/778.  =  c^AT 

The  driver  gas  is  initially  at  room  temperature  and  it  is  assumed  that  the 
small  temperature  decrease  during  the  cycle  is  not  influenced  by  heat  transfer. 

Piston  motion  is  evaluated  by  applying  a  force  balance  across  the 
piston,  taking  into  account  the  frictional  drag. 


F  =  MpD  '  Pt)  “  D 

where  pD  and  p-p  are  the  respective  pressures  of  the  driver  and  test  gases,  and 
D  is  the  frictional  drag  of  the  piston  which  is  expressed  in  terms  of  piston 
velocity  by: 


Piston  acceleration,  velocity  and  travel  follow: 


a 


P 


_F 

m 


AV  =  a^At 
P  P 


Axp  =  -^At2  +  v0At 


Test  Gas  Compression 

The  compression  of  the  test  gas  occurs  as  a  result  of  piston 
motion.  Energy  that  is  added  by  virtue  of  the  compression  is  calculated 
from  the  conservation  of  energy  equation.  The  work  done  by  the  piston  on  the 
test  gas  during  this  compression  is  one  term  in  this  conservation  of  energy 
equation.  Other  terms  include  heat  loss  to  the  wall  of  the  tube,  which  becomes 
important  as  the  gas  temperature  rises.  The  other  important  equation  is  con¬ 
servation  of  mass.  The  test  chamber  is  not  a  closed  chamber  but  contains 
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an  exhaust  port  at  the  downstream  end  where  the  test  sample  is  located.  Thus, 
test  gas  is  allowed  to  flow  from  the  test  chamber  through  the  test  sample. 
Therefore,  the  mass  in  the  system  is  not  constant  by  virtue  of  mass  and 
energy  flow  from  the  driven  tube  and  plenum  chamber  through  the  test  sample. 
Terms  in  the  conservation  of  energy  and  mass  equations  reflect  this  mass 
and  energy  loss. 

The  equation  of  state  that  applies  to  the  test  gas  is  the  Van  der 
Wall  equation  which  includes  terms  to  express  the  nonlinear  relationship 
between  pressure  density  and  temperature.  The  terms  for  this  equation  are 
determined,  as  mentioned  previously,  by  a  linear  averaging  of  the  mole 
fraction  of  the  test  gas  constituents. 

The  test  gas  specific  heat  is  also  assumed  to  be  for  non-perfect 
gas  and  is  expressed  in  terms  of  a  linear  function  of  the  gas  temperature. 

The  coefficients  in  this  expression  are  also  linear  averages  of  the  mole 
fraction  of  the  test  gas  compositions.  The  test  gas  specific  heat  is  expressed 
in  terms  of  a  secant  function  in  which  the  product  of  the  specific  heat  and 
the  temperature  yield  the  internal  energy.  This  is  contrary  to  a  normal 
expression  of  specific  heat  whereby  it  is  a  tangent  function  so  that  the 
integral  of  the  product  of  specific  heat  and  temperature  yield  the  internal 
energy.  The  technique  used  here  provides  rather  simple  yet  effective  means 
for  evaluating  the  internal  energy  in  a  finite  difference  scheme  with 
many  time  steps. 

The  Van  der  Waal  equation  of  state  used  for  the  test  gas  is 


where 


ai 

ei' 


27 


RT 


ci 


ai 

fTT  =  27pci’ 


v  is  the  specific  volume,  and  Tci  and  pci  are  the  critical  temperature  and  pres¬ 
sure  for  the  ith  gas  constituent,  a  and  £3  are  the  average  quantities  based  on 
the  mole  fraction  of  the  ith  constituent. 


The  specific  heat  at  constant  volume  is  defined  by 


c 


v 


f (c  +  E.X  .  c  T. (T  -  460)ei)T  -  ^ 

V  V  VO  1  gl  V  1 


778; 


1_ 

Ti 


where  cvo  and  CyTj  are  the  intercept  and  slope  of  the  temperature-dependent 
specific  heat,  Xgi  is  the  mole  fraction  of  the  ith  constituent,  and  ei  is  the 
temperature  exponent  in  the  curve  fit  relationship. 

The  test  chamber  mass  balance  is  given  by 

ntj.  =  mTo  -  Am 
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where  Am  is  the  mass  flow  through  the  test  sample.  Initially  this  flow  is 
assumed  to  be  negligible  and  the  test  chamber  and  barrel  volumes  are  lumped 
together.  When  the  projectile  velocity  exceeds  100  ft/sec,  the  calculation 
of  flow  through  the  test  sample  is  initiated.  This  computation  involves 
determination  whether  sonic  or  subsonic  flow  conditions  exist  within  the 
test  sample.  The  sonic  static  pressure  is  given  by: 


p* 


Px  (1  + 


Y  -  1,-W  -  1) 
2  ’ 


where  p^  is  the  test  chamber  pressure. 

If  p*  is  greater  than  the  barrel  pressure,  P2,  (downstream  from  the  test 
sample)  then  sonic  conditions  exist  and 

Am  =  p*A12(^f)1/2At, 

where  Aj2  is  the  test  sample  flow  area,  T*  is  the  sonic  static  temperature, 
and  At  is  the  computation  time  interval. 

If  subsonic  conditions  exist. 


Am  =  8.02  A12((P1  -  p2)  ^)1/2At, 

which  is  the  equation  for  flow  through  a  venturi  in  terms  of  the  upstream 
and  downstream  line  pressures. 

The  change  in  internal  energy  in  the  chamber  over  the  calculation 
time  interval  is  given  by 

PiVxp 

AE1  =  —778  P  -  CVAmTlY 

which  represents  the  compression  work  due  to  the  piston  and  the  loss  in 
enthalpy  due  to  flow  through  the  test  sample.  Heat  transfer  losses  in  the 
driven  tube  and  chamber  are  not  included  in  this  analysis  at  present .  The 
gas  chamber  temperature  is  defined  by 


where  is  the  current  value  of  internal  energy. 
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The  test  sample  is  a  straight  channel  with  a  rounded  entrance. 

The  flow  through  the  channel  is  computed  by  either  sonic  or  subconic  condi¬ 
tions  depending  on  the  pressure  at  the  inlet  and  outlet  to  the  sample.  The 
calculations  provide  the  static  flow  conditions  of  pressure,  temperature,  and 
density  in  addition  to  the  flow  velocity  over  the  surface  of  the  sample.  These 
conditions  are  in  turn  used  in  the  equation  that  expresses  turbulent  heat 
flux  to  a  flat  plate.  The  heat  flux  is  computed  and  then  summed  to  yield  a 
current  level  of  total  heat  input.  The  heat  flux  to  the  surface  is  also 
applied  to  an  unsteady  heat  conduction  routine  by  which  the  surface  tenpera- 
ture  and  the  temperature  distribution  in  the  test  sample  are  evaluated. 

These  calculations  are  all  performed  within  the  same  time  step  of  the  overall 
finite  difference  calculation. 

The  technique  used  to  calculate  heat  flux  to  the  test  sample 
surface  requires  the  flow  velocity,  density,  and  viscosity.  The  density  and 
viscosity  are  evaluated  on  the  basis  of  a  reference  temperature  to  take 
the  temperature  profile  resulting  from  the  boundary  layer  velocity  distribution 
and  sample  surface  temperature  into  account. 

Mach  number  of  flow  through  the  test  sample: 


where  p  is  the  sonic  static  pressure  if  M  s  1  or  the  downstream  pressure,  p2 
if  the  flow  is  subsonic. 

Free  stream  static  temperature: 

Te  =  Tx(l  +(^~j~)M2)-1 

Free  stream  velocity: 

Ue  =  [2gJ  YCV(T1  -  Te)]1/2 


Reference  temperature: 


Tref 


<Te  +  W/2 


0.11  Ue2 

gJ  CVY 


SAMP 


where  T, 


is  the  test  sample  surface  temperature. 


Viscosity: 


=  7.0  x  10-7  T 


1.5 


ref 


<Tref  + 


198) 


-1 


Density: 

pref  ■  PS  <6PS  *  RTrafJ 


Turbulent  flat  plate  heat  flux  at  the  location  of  the  in-wall  thermocouple: 


Q  ■  0.052(PrefUe)°'V'2  CVT  (T, 


T  ) 
SAMP; 


The  sample  surface  temperature  is  determined  from  the  one-dimensional 
unsteady  state  heat  conduction  equation. 


9T  _  a^r 

at  aa  2 

X 


with  the  surface  boundary  condition. 


for  X  =  0,  where  Q  is  the  heat  flux  to  the  sample  surface.  This  equation  does 
not  consider  the  effects  of  cylindrical  geometry. 

A  finite  difference  technique  using  a  geometrical  node  grid  spacing 
was  incorporated  into  the  Calspan  code  to  solve  the  unsteady  heat  conduction 
equation.  The  general  finite  difference  relationship  is  given  by 


2a _ 

(F+1)Ax2._j 


(T 


i-1 


+  —  T 
F  i+1 
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where  a  is  the  thermal  diffusivity 
T  is  the  temperature  rise 
t  is  the  time  interval 

Ax^_i  is  the  thickness  of  the  i-lst  grid,  and 

F  is  the  geometrical  multiplier  with  the  thickness  of  the  ith 
grid  being  F  times  that  of  the  i-lst  grid. 

The  exposed  surface  boundary  condition  is  satisfied  by  first  determining  a 
fictitious  temperature  in  free  space, 
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2q 


+  T 
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This  temperature  is  then  used  to  establish  the  surface  temperature  rise  by 


5 

At 


ST2  <To  -  2T1  +  V 


Barrel  Flow 


The  flow  through  the  barrel  is  expressed  as  an  input  of  mass  and 
energy  to  the  volume  between  the  test  sample  and  the  projectile.  As  mass 
and  energy  are  accumulated,  this  is  expressed  in  terms  of  pressure  and 
temperature,  which  in  turn  provides  the  accelerating  force  for  the  projectile. 

In  this  calculation,  the  quasi-steady  assumption  of  the  previous  calculation 
is  relaxed  by  allowing  pressure  acting  on  the  base  of  the  projectile  to  be 
modified  according  to  the  Mach  number  of  the  flow  at  the  base  of  the  projectile. 
In  this  way,  the  unsteady  expansion  effects  of  the  flow  through  the  barrel  is 
taken  into  account.  The  equation  of  state  and  the  basic  energy  and  mass 
conservation  equations  are  the  same  as  for  the  test  gas  in  the  driven  tube. 

These  equations  are  used  to  define  "2"  conditions  in  the  barrel. 


Projectile  motion  is  calculated  by  an  approximate  technique  that 
involves  determination  of  the  flow  Mach  number  at  the  projectile  base. 


Specific  enthalpy: 


2  Y 


M 


where  Mg  is  the  mass  of  gas  contained  in  the  barrel. 


Static  enthalpy  at  the  projectile  base: 
d  V  2 

H2  =  2.5  x  10  H2  -  -|- 


where  V 

P 


is  the  projectile  velocity 
Mach  number: 


M  =  V  {(y  -  1)  H  }"1/2 
‘  00 


The  assumption  is  made  that  the  pressure,  p2>  is  not  the  total  pressure,  but 
is  a  static  pressure  at  the  barrel  origin. 
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It  is  further  assumed  that  the  difference  in  pressure  between  this  location 
and  the  projectile  base  is  equal  to  the  difference  between  the  pressure  and 
the  total  pressure.  This  is  expressed  as  a  Mach  number  function: 


P  -  Ap  Yj_  9  -Y/Y-l 

ynrv  -  a  +  2  mj  . 


This  allows  the  static  pressure  at  the  projectile  base  to  be  expressed  in 
terms  of  the  barrel  origin  pressure, 

From  this  the  projectile  acceleration: 


a 

pr 


32.2  (P2  -  p  )/W 

oo  * 


where  p^  is  the  projectile  resistance  to  motion,  and  W  is  the  projectile 
weight,  with  Pr  =  Pr  +  Pr  x  P 

o  x  pr‘ 

Projectile  velocity  change: 

A V  «  a  At 

pr  pr 


Projectile  displacement: 

At2 

AX  =  V  At  +  a 
pr  pr  pr  2 
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c 

c 

c 

C  LIST  OF  CONSTANTS  AND  VARIABLES  USED  IN  THIS  SHOCK  TUBE  GUN  PROGRAM 
C 

C  DRIVER  GAS  PARAMETERS 

C  PCH  INITIAL  DRIVER  PRESSURE  PSIA 

C  VOLO  INITIAL  VOLUME  OF  PROPELLING  GAS  -  FT**3 

C  TO  INITIAL  PROPELLING  GAS  TEMPERATURE  -  DEG  R 

C  CVDVR  DRIVER  GAS  SPECIFIC  HEAT  -  BTU/LBM-DEG  R 

C  RDVR  DRIVER  GAS  CONSTANT  -  FT-LBF /LBM-DEC  R 

C 

C  PISTON  PARAMETERS 

C  AT  PISTON  AREA  -  FT**2 

C  KP1  DRAG  RATE  CONSTANT  -  LBF/FT 

C  VELP  INITIAL  PISTON  VELOCITY  -  FT/SEC 

C  WSHELL  PISTON  WEIGHT  -  LBM 

C  SHELRV  PISTON  MASS  DURING  RETURN  -  LBM 

C  XPMAX  MAX  PISTON  DISPLACEMENT  -  FT. 

C 

C  TEST  GAS  PARAMETERS 

C  P20  INITIAL  TEST  GAS  PRESSURE  -  PSIA 

C  T20  INITIAL  TEST  GAS  TEMPERATURE  -  DEG  R 

C  R2  TEST  GAS  CONSTANT  -  FT-L BF/L BM-DEG  R 

C  CV2  SPECIFIC  HEAT  OF  TEST  GAS-  BTU/LBM-DEG  R 

C  CVT 2  TEST  GAS  SPECIFIC  HEAT 

C  CVEXP  CURVE-FITTING  EXPONENT 

C  AI.FA2  VANDERWAAL  CONSTANT 

C  BETA2  VANDERWAAL  CONSTANT 

C  GAM2  TEST  GAS  RATIO  OF  SPECIFIC  HEATS 

C 

C  TEST  VOLUME  PARAMETERS 

C  A 2  TEST  VOLUME  AREA  -  FT**2 

C  D3  DIAMETER  OF  TEST  VOLUME  EXIT  -  IN 

C  CD2  FLOW  COEFF.  OF  TEST  VOLUME  EXIT 

C  VOL20  INITIAL  VOLUME  OF  TEST  GAS  IN  TUBE  -  FT**3 

C  VOL2F  FINAL  TEST  VOLUME  -  FT**3 

C 

C  PROJECTILE  PARAMETERS 
C  BARL  CARREL  LENGTH  -  FT 

C  BORED  ROBE  DIAMETER-FT 

C  PSTART  PRESSURE  ABOVE  WHICH  PROJECTILE  IS  ALLOWED  TO  MOVE  -  LB/IN**2 

C  RESO  Y-INTERCEPT  OF  PROJECTILE  RESISTANCE  FUNCTION  -  LBM 

C  REGS  SLOPE  OF  PROJECTILE  RESISTANCE  FUNCTION  -  LBM/FT 

C  RESC  CONSTANT  RESISTANCE  TERM  AFTER  FUNCTION  L.T.  RESC  -  L3M 

C  VPPvOJ  INITIAL  PROJECTILE  VELOCITY  -  FT/SEC 

C  WPROJ  PROJECTILE  WEIGHT  -  LBM. 

C  XPROJ  INITIAL  PROJECTILE  DISPLACEMENT  -  FT 

C 

C  TIME  PARAMETERS 
C  DELT  TIME  INCREMENT  -  SEC 

C  PRCI  NUMBER  OF  TIMES  BETWEEN  PRINT  INTERVALS 

C  TFO  LIMITING  TIME  -  SEC 

C 
C 
C 

READ  <5,101)  NUMBER 
DO  99  NM-1, NUMBER 

DIMENSION  RG<9),CVG<9),CVTG(9),CVEXPG<9),ALFAG<9),BETAG<9),XG<9) 
DIMENSION  G1  ( 9  )  ,G2< 9  )  ,G3( 9  )  ,G4< 9 >  ,G5< 9 > 

DIMENSION  TNEWt  40 ) ,TSAMP<  40 ) , DELX( 40  > ,D!<NEW<  40  )  ,DTEMP<  40  > 
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REALM  M2,M20,M1  ,KP1 

DATA  VOLO/31 .  5/ .CVDVR/0 . 1 76/ , RDVR/55 . 2/ 

DATA  AT/O .3068/, WS HELL/150. 0/.SHELRV/150.0/.XP MAX/81.0/ 

DATA  A2/0.3068/ , D3/0 . 5/ ,CD2/0 . 75/ , V0L2F /O . 00223/ 

DATA  BARL/ 1 5 . 0/ , BORED/O . 0984/ , RES0/7200 . / , RESS/4320 . / , RESC/7200 . / 
DATA  WPROJ/O . 25/ i VOL 1F/0.082/.DOR12/0.5/ 

JJ-1 

KK-1 

NN-0 

OFLUX-O.O 

UE-0.0 

R2  =  0 . 0 

CV2-0.0 

CVT2-0.0 

CVEXP-0.0 

ALFA2-0.0 

BETA2-0.0 

GAM2-0.0 

P20-14.7 

T20-530.0 

XPROO-O.O 

VELP-0.0 

VPROO-O.O 

VOL20-0.0 

D0R12-0.5 

100  FORMAT  <  8F 1 0 . S  ) 

101  FORMAT  <  II  ) 

102  FORMAT  (  F5 . 2 , 5A3 , 6E 1 0 . 3  ) 

103  FORMAT  (12H1  INPUT  DATA.T48, ' TEST  N0.',T59,I1> 

104  FORMAT  < ////T6 . ' TEST  GAS  MIXTURE  BY  MOLE  FRACTION'/) 

105  FORMAT  ( /T6 . F6 . 3 , T1 7 . 5A3 > 

READ  (5,100)  PCH,T0,KP1 , P START , DELT , PRC  I ,TFO , FACTOR 
READ  (5,101)  NG 
WRITE  (6.103)  NM 
WRITE  (6.104) 

DO  3  ID-1 ,NG 

REAO  (5,102)  XG( !D),G1( I D ) , G2 ( I D ) , G3( ID),G4( ID),G5( ID > ,RG( ID > ,CVG( 
♦  ID), CVTG(  ID ) , CVEXPG( ID ) , ALFAGf ID), BETAG( I D  ) 

WRITE  (6,105)  XG( I D  ) , G 1 ( I D  )  , G2( I D  >  ,  G3(  I D ) , G4( ID),G5( ID) 

R2  -  R2  +  XG( ID )*RG( ID) 

CV2  -  CV2  ♦  XG( ID )*CVG( ID  ) 

ALFA2  -  ALFA2  ♦  XG( ID )*ALFAG( ID  ) 

BETA2  -  BETA2  ♦  XG( I D  )*BETAG( ID  > 

3  CONTINUE 

GAM2  -  1.  +  ( R2/(  778 . *CV2  > ) 

ABORE  -  0.785*  BORED**2 

CV-CVDVR 

DELTO-DELT 

DOR  1 2-DOR 1 2/ 1 2 . 0 

A12»0.7854*DOR12**2 

IPRC-IF1X(PRCI  ) 

I  PR  I  NT-0 
PCH-PCH* 1 44 . 

PCHP-PCH/ 1 44 . 

P20-P20*144.0 
PARTIM  -  0. 

PSTART-PSTART*144. 

QTOT-O.O 

RGAS-RDVR 

CP-RGAS/778.+CV 

GAM-CP/CV 
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GM  -  PCH/<  RGAS*TO )*VOLO 
E-GM*( CP-RGAS/778 .  )*T0 
TIME  -  0. 

VMAX-XPMAX*AT 
XP -XPMAX 

C  INITIAL  CONDITIONS  -  UPSTREAM  -  FROM  PISTON  FACE  TO  TEST  SAMPLE 
P1-P20 

VOL 1-VOL 1F+XP*A2 
T 1 -T20 

M1»P1*V0L1/R2/T1 

CVT2-0.0 

DO  1  0-1, NG 

CVT2-CVT2*XG< J  >*CVTG(J >*< T1 -460 . 0 >**CVEXPG< J  ) 

1  CONTINUE 

E 1-MI "I CV2+CVT2 )*T1 
RHOl-Ml/VOLl 
VV1-1 .0/RH01 

C  INITIAL  CONDITIONS  -  DOWNSTREAM  -  FROM  TEST  SAMPLE  TO  PROOECTILE  BASE 
P2-P20 

VOL2»VOL2F+XPROJ*ABORE 

T2-T20 

M2-P2*VOL2/R2/T2 

CVT2-0.0 

DO  2  J-l , NG 

CVT2-CVT2+XG( J)*CVTG(J)*(T2-460.0  )**CVEXPG( 0  ) 

2  CONTINUE 
E2-M2*(CV2+CVT2)*T2 
RH02-M2/V0L2 
VV2-1 ./RH02 

WRITE  (6,110) 

110  FORMAT  < ///T6, ' INITIAL  CONDITIONS  FOR  TEST') 

WRITE<  6, 120  )P20,PCH,AT,T20,TO,KP1 , R2 , RDVR , WSHELL , C V2 , CVDVR , SHELRV , 
+CVT2. VOLO, XPMAX, CVEXP, PSTART.ALFA2 , BARL , BETA2 , A2 , BORED, GAM2.D3, 

+  DOR12,CD2,RESO,DELT,VOL20,RESS,PRCI, VOL2F , RESC , TFO , VOL  IF , WPROJ 
120  FORMAT  <//T6,‘P20  -  1 ,T1 7 , F 1 3 , 6 ,T3 1 , • PSF '  ,T48,,PCH  - ' ,T59 , F 13 .6 , 
♦T73 , 1 PSF '  , T90 , 'AT  -  1 , T1 0 1 , F 1 3 . 6 ,T1 15 , ' FT**2 ' / /T6 , ' T20  «‘,T17, 
+F13.6.T31, 1 DEGR ' ,T48, 'TO  - ' , T59 , F 1 3 . 6 , T73 , ' DEGR ' ,T90, 'KP1  •* .T101, 
+F13.6.T115, 1 LBF /FT  1 //T6 , 'R2  -  1 ,T 1 7 . F 1 3 . 6 , T3 1 , '  FT-LBF /LBM-DEGR 1 , 
♦T48, 'RDVR  , T59, FI 3. 6, T73, ’FT-LBF /LBM-DEGR' ,T90, 'WSHELL  -' .T101. 
+F13.6.T115, ' LBM ' //T6 , ' CV2  - ' , T1 7 . F 1 3 . 6 . T3 1 , ' BTU/LBM-DEGR ’ , T48 , 
♦'CVDVR  .T59.F13. 6, T73, 'BTU/LBM-DEGR' ,T90, 'SHELRV  - ' , T1 0 1 . F 1 3 . 6 , 
+T115, 'LBM'//T6, ' CVT2  - ' , T 1 7 , F 1 3 . 6 , T48 ,  ' VOLO  -  * , T59 , F 1 3 . 6 ,T73 . 
+'FT**3' ,T90, 'XPMAX  - ' , T10 1 ,F 1 3 . 6 ,T1 15 , ' FT ' //T6 , ' CVEXP  -',T17, 
+F13.6.T48, 'PSTART  - ' , T59 , F 1 3 . 6 , T73 , ' PSF '  //T6 , ' ALFA2  - ' , T1 7 , F 13 . 6 , 
+T90 , ' BARL  ■' ,T101,F13.6,TI15, ' FT ’ //T6 , ' BETA2  • ' ,T17 , F 13 .6 ,T48, 

+  'A2  -' ,T59,F13.6,T90, 'BORED  - ' , T 1 0 1 , F 1 3 . 6 , T1 1 5 , ' FT ' //T6 , ' GAM2  »', 

+  T17.F13.6.T48, ' D3  «' , T59 , F 1 3 . 6 , T90 , • DOR  1 2  -',T101,F13.6,T115,'FT'/ 
+/T48 , ' CD2  -' ,T59,F13.6,T90, 'RESO  - ’ ,TI 01 , F 1 3 , 6 , T1 15 , ' LBM ' //T6 . 

+ 'DELT  *' ,T17,F13.6,T31 , 'SEC' ,T4B, 'VOL20  , T59 , F 1 3 . 6 , T73 , 'FT**3' , 
+T90, 'RESS  *' ,T101,F13.6,T115, ' LBM/FT ' //T6 , 'PRC I  - ' , T1 7 , F 1 3 . 6 , T48 , 

♦ ' VOL2F  -' ,T59,F13.6,T73, 'FT**3' ,T90, 'RESC  = ' , T1 0 1 , F 13 . 6 ,T1 1 5 , ' LBM ' 
+  //T6 , 1 TFO  ■' ,T17,F13.6,T31, 'SEC' ,T48, ' VOL  IF  - ' , T59 , F 1 3 . 6 ,  T73 , 

+  ' FT--3 ' , T90 , ' WPROO  - ' , T1 0 1 , F 1 3 . 6 , T1 1 5 , ' LBM ' ///// ) 

DEPTH  -  0.5/12.0 
XK  -  19.3/3600. 

ALPHA  »  XK/57.6 
F-l  .3 

DELX0-SQRT( ALPHA*DELT/0 . 25  ) 

SUMX-DELXO 
DELX< 1  )-DELX0 
DO  6  1-1,40 
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TSAMP ( I  )-T0 
TNEW< I >-T0 

6  CONTINUE 

OO  7  1-1,39 

DELX< 1*1 >-DELX( I >*F 

IF< SUMX.GE .DEPTH >  GO  TO  7 

SUMX-SUMX+DELX( I > 

KTEMP  -  I 

7  CONTINUE 

CONST 1 -ALP HA* 2 . 0/( l . 0*F  > 

CONST2-I 1 .  0  +  F  )/F 
CONST3- 1 . 0/F 

TZIP»QFLUX*2 . 0*DELX( 1 )/XK+TSAMP<  2  ) 

NAME L I  ST/FIRST/DE LX, KTEMP, TZ IP 
WRITE  (6, FIRST) 

10  CONTINUE 
T-TO 

RGAS-S7 .4*.333E-5*PCH 
CV-0. 175*. 183E-4*T 
CP-RGAS/778 . +CV 
GAM-CP/CV 

IF <  VELP . LT. -0 . 1 >  WSHELL-SHELRV 

C  COMPUTE  THE  DYNAMICS  AND  THERMODYNAMICS  OF  PISTON  MOTION. 
24  XD-XPMAX-XP 

DRAG-KP 1*ABS( VELP > 

FOVM-<  (  PCH*AT-P 1 *A2 (-DRAG )/WSHELL*32. 17 
26  DVELP  -  FOVM  *  DELT 
28  VELP  -VELP  +DVELP 

IF IXP.LE. 0.00001 . AND. VELP. GT. 0.0)  VELP-0.0 
DXP-DELT*< VELP+F0VM*DELT/2. > 

XP-XP-DXP 

I F ( XP . LE .  .00001  >  GO  TO  99 
WORK  -  PCH  *  AT  *  DXP/778 . 

E  -  E  -  WORK 

30  VOL=VOLO*AT*XD 
RHOCH  -  GM/VOL 

T  =  E/<  < CP-RGAS/778  .  >*GM) 

PCH  -  RHOCH*RGAS*T 
C  COMPUTE  TEST  GAS  PROPERTIES. 

CVTX-0.0 
CVTY-0.0 
DO  31  J-l.NG 

CVTX-CVTX+XGI 0 )*CVTG( 0 )*(T2-460.0  )**CVEXPG( J  ) 
CVTY-CVTY  +  XGI3 )*CVTG< J )*< T1 -460 . 0 >**CVEXPG< J  ) 

31  CONTINUE 

CVY  =  ( <  CV2+CVTY  )*T 1 -ALFA2*RH01 /77B . 0  >/Tl 
CVX-< (CV2+CVTX  )*T2-ALFA2*RH02/778 . 0  )/T2 
RESIST  -  RESO  +  RESS  *  XPROJ 

32  IF< RESIST. LT.RESC  )  RESIST-RESC 
H2-E2*GAM2/M2 

H2INF-25000. *H2-VPROO**2/2 . 

XM2 I NF -VPROO/SQRTl <  GAM2- 1 .0)*H2INF  ) 

FM-(  1 .0+IGAM2-1 ,0)/2.0*XM2INF**2 >**< -GAM2/I GAM2-1 . 0  ) ) 
DELP»P2*( 1 .0-FM)/< 1 .0*FM> 

33  P20-P2+DELP 
P2INF-P2-DELP 

PFORCE  ■  ABORE  *  ( P2 INF-RES  I  ST > 

I F ( PFORCE  .  LT .  0 . 0  )  PFORCE-O.O 

IFIP2.LT. PSTART  )  PFORCE-O.O 

APROO  =  P FORCE /WPROJ *32 . 1 7 

DXPROO  -  VPROO  *  DELT  ♦  APROO’  *  DELT**2/2. 
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35  XPROO  -  XPROO  +  DXPROJ 

I F < XPROJ . GE . BARL >  GO  TO  99 
v  .JO  ■  VPROJ  +  APROJ"  DELT 
VOL  1 =VOL 1F+XP*A2 
VOL2«VOL2F+XPROJ*ABORE 

36  P1AVE-P1 
XM1-M1 
XM2=M2 
IPASS-0 
P 1 SAV«P 1 
P2SAV-P2 

37  CONTINUE 
IP/'SS-IPASS  +  l 

P STAR«P 1 *< l .0+<GAM2-l .0  1/2.0  )•«< -GAM2/I GAM2- 1 .0)  > 

TTOT-T1 

PS=*PSTAR 

PTOT-P1 

I  F( P20 .GT . P 1 >  GO  TO  47 
I F ( PSTAR . GE . P20  )  GO  TO  38 
PS»P20 

DELTM1”8.02*A12*SORT(<P1-P20)*P1/R2/T1  )*DELT 
GO  TO  49 

38  TSTAR»Tl*2.0/( GAM2+ 1.0) 

DELTM1*PSTAR*A1 2*SQRT( GAM2*32 . 2/R2/TSTAR  >*DELT 
GO  TO  49 

47  CONTINUE 
PSTAR-PSTAR*P20/P1 
TTOT»T2 
PS-PSTAR 
PTOT»P20 

IF<  PSTAR. GE. PI  )  GO  TO  48 
PS  =  P1 

UELTM1=-8.02*A12*SQRT( ( P20-P1  )*P20/R2/T2  >*DELT 
GO  TO  49 

48  TSTAR«T2*2.0/<  GAM2+ 1 . 0  > 

DELTM1 *-PSTAR*AI 2*SQRT( GAM2*32 .2/RZ/TSTAR)*0£LT 

49  CONTINUE 
XM1-M1-DELTM1 
XM2-M2+DELTM1 
RH01-XM1/V0L1 
RH02-XM2/V0L2 
VV1-1 .0/RH01 
VV2-I ./RH02 

WORK»-ABORE*DXPROO*P2 I NF/778.0-P2/Z88. “XPROO *3 . 1416“BORED“DELT 
C  P2/288 .  { P2  IS  IN  PSF >  IS  ASSUMED  APPROXIMATELY  EQUAL  TO  HEAT  FLUX. 
I F  <  VPROO . LT . 1 00 .  )  GO  TO  58 
DELH=CVY*DELTM1 *T1 “GAM2 

50  I F ( 0ELTM1 . LT . 0 . 0  )  DELH-  CVX*DELTM1“T2*GAM2 
EX=E2+WORK+OELH 

T2-EX/CVX/XM2 

PX  =  <  R2*T2 )/( VV2-BETA2 )+ALFA2/VV2**2 
P  2 “ < PX  +  P2SAV  >/2 . 0 
EV-E1-DELH“P1AVE/77S.0*A2*DXP 
TI aE V/CVY/XM1 

PV-(R2*T1  >/< VV1-BETA2)+ALFA2/VV1“*2 
Pl-(PY+PISAV>/2.0 

54  DELP1-PY-P1 
P1AVE-PH-DELP1/2.0 
P20-P2+0ELP 

55  IF( IPASS.LT. 3)  GO  TO  37 
El-EY 
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E2-EX 
Pl-PY 
P2-PX 
Ml “XM1 
M2  ■  XM2 

XMOR-SQRK  ( (  PTOT/PS  >**<  <  GAMZ-1 . 0 I/GAM2  >- 1 . 0  >»2 . 0/< GAM2- 1 . 0  )  > 
TE»TTOT/< 1 .0  +  <GAM2-l . 0 1/2 . 0*XMOR**2 > 

U£  *SQRT <  S . OE  4*CVY*GAM2*( TTOT-TE > ) 

TREF “(  TE+TSAMPl  1 > > /2 . 0»4 . 4E-6*UE**2/CVY/GAM2 
VIS»7.0E-7*TREF**1.5/<TREF+198.0) 

RHOREF“PS/< BETA2*PS+R2*TREF  ) 

QFLUX“0.052*<  RHOREF  *UE  ) **0 . 8*VI S**0 . 2*CVY*GAM2*< TTOT-TSAMP ( 1 ) ) 

I F ( TSAMP (  1 I.GT.3100.0  >  QFLUX“QFLUX*< TTOT-31 00 . 0 >/< TTOT-TSAMP ( 1 ) ) 
QFLUX  *  QFLUX*FACTOR 
TZIP“OFLUX*2.0*DELX( 1 I/XK  +  TSAMP <  2 > 

OTEMP  < 1 >-ALPHA/(DELX< 1 >**2 >*( TZIP-2 . 0*TSAMP< 1 )*TSAMP(2>  > 

DO  56  K-2.KTEMP 

DTEMP<K>-C0NST1/(DELX<K-1 >**2  >*< TSAMP <K-1 >-C0NST2*TSAMP< K >+ 

*  C0NST3*TSAMP(K+1 > ) 

56  CONTINUE 

DO  57  K“ 1 , KTEMP 

TSAMP ( K  > *TSAMP ( K )+DTEMP<  K  )*DELT 

57  CONTINUE 

TSAMPl KTEMP  +  1 )-TSAMP<  KTEMP-1  ) 

OTOT-QTOT+QFLUX*DELT 
I F {  P20.GT.P1  >  UE  —  UE 
GO  TO  60 

58  CONTINUE 

VOL =VOL 1 +VOL2 
XMTOT =M1 +M2 

ESUM-El+E2+Pl»A2*DXP/778.0  +WORK 

Ml =XMTOT*VOL 1 /VOL 

M2 “XMTOT* VOL 2/ VOL 

El“ESUM*VOLl/VOL 

E2“ESUM*VOL2/VOL 

T1-ESUM/CVX/XMT0T 

T2-T1 

VVT-VOL/XMTOT 

P1-(R2*T1 >/( VVT-BETA2  >*ALFA2/VVT**2 
P2-P1 

RH01-M1/V0L1 
RH02-M2/V0L2 
VV1-1 .0/RH01 
VV2-1 .0/RH02 
60  CONTINUE 

IF  ( VOL  1 .GT . 0 . 2*VMAX  )  GO  TO  80 
I F ( VOL  1 .LE.0.2*VMAX)DELT-0.1*DELTO 
IF(VOL1.LE.O. 2* VMAX .AND.KK.EQ. 1 >  GO  TO  62 
IF< VOL1 .LE.0.02*VMAX)DELT“0.01*DELTO 
IF(VOL1.LE.0.02*VMAX. AND.KK.EQ. 2)  GO  TO  62 
GO  TO  80 

C  REVISE  GRID  SIZE  DUE  TO  CHANGE  IN  TIME  STEP. 

62  0ELX0-SQRT(ALPHA*DELT/0.25> 

SUMX-DELXO 

DXNEW( 1  I-DELXO 

DO  64  1-2,40 

DXNEW( I  )“DXNEW( 1-1 >*F 

I F < SUMX.GE .DEPTH  )  GO  TO  64 

SUMX*SUMX  +  DXNEW< I  ) 

KNMAX  -  I 
64  CONTINUE 


77 


XKOLD-DELX!  1  ) 

XKNEW-0.0 

K-l 

KNEW* 1 

A  =  !TZIP-TSAMP<2  >  >/<  2 . 0*DELX!  K ) > 

B-(TSAMP< 1 )-TSAMP !2)-A*DELX!K))/DELX!K)*'*2 
TN£W< 1  )-TSAMP( 1  ) 

66  KNEW-KNEW+1 
XKNEW-XKNEW+DXNEW! KNEW-1  ) 

TNEW!KNEW)=TSAMP! 1 >-! A*XKNEW+B*XKNEW**2 > 

Ir! XKNEW.LT. XKOLD  )GO  TO  66 

67  K-K+l 

XKOLD«XKOLD+DELX( K ) 

A»( TSAMP<  K- 1 >-TSAMPI K+l  )  >/! ( F  +  l . 0  >«DELX! K- 1 > ) 

B  =  !  TSAMP  <  K  )-TSAMP  !  K  + 1  )-A*DELX(K))/DELX!K)**2 

68  KNEW-KNEW+1 
XKNEW-XKNEW+DXNEW! KNEW-1 ) 

XDEL-XKNEW-XKOLD+DELX!  K ) 

TNEW( KNEW ) ’TSAMP <K)-!A*XD£L+B*XDEL**2  ) 

I F ( XKNEW. LT . XKOLD .AND .  KNEW+ 1 . LE . KNMAX >  GO  TO  68 
I F<  K+l . LE  . KTEMP . AND . KNEW+1 . LE . KNMAX  >  GO  TO  67 
DO  70  1=1 .KNMAX 
DELX< I )=DXNEW( I > 

TSAMP ( I >-TNEW( 1  ) 

TNEW< I  )=T0 
70  CONTINUE 
KTEMP-KNMAX 
KK  =  KK  + 1 

NAMEL I  ST/SECOND/DEL  X. KTEMP, T2IP, A, B 
WRITE  <6. SECOND) 

80  CONTINUE 

C  PRINT  OUT  COMPUTED  RESULTS. 

IF  ( 00. EQ. 1)  GO  TO  88 
87  NN-1 

PARTIM-PARTIM+DELT 
I F<  T IME .GT . TFO  )  GO  TO  99 
I F (  ^RINT.LT. IPRC  >  GO  TO  90 
GO  TO  89 

187  FORMAT^  f isHl ^TABULATED  OUTPUT) 

WRITE  <6, 188) 

188  FORMAT!  /  //T6, 'TIME' ,TZ0, 'PI' ,  T34 , ‘ T1 ' ,T4B , 'Ml 1 ,T62, 'VOLl 1 ,T76, ’VEL 
♦P' , T90 , ' VPROJ ' , T1 04 , 'WORK' ,T118, 'QFLUX'/T6, ' PCH' ,T20, ' P2 ' , T34 , ' T2 1 
+  .T48, 'M2' ,T62, 1  VOLZ  1 ,T76, 'XP' ,T90, ' XPROJ ' ,T104 , ' UE ' ,T1 18, ' QTOT 1 / 
♦T4, 'TSAMP! 1  )  ' ,T18, 'TSAMP! 2) ' , T32 , ' TSAMP! 3 ) ' ,T46 , ' TSAMP! 4 ) * , 

♦T60 , ' TSAMP! 5  ) ' . T74 . ' TSAMP ! 6  > ' , T88 , ' TSAMP! 7 ) ' , T1 02 , ' TSAMP ! 8 ) ' . 

♦  T1 16 , ' TSAMP! 9  ) ' /  /  ) 

00=00+1 

89  PCHP-PCH/144. 

P1P=P 1/144.0 
P2P-P2/144. 

TIME-TIME+PARTIM 
PARTIM  ■  0. 

WRITE! 6, 190)  T IME . P 1 P , T1 ,Ml , VOL  1 ,VELP , VPROO .WORK .QFLUX ,PCHP , P2P ,T2 
+  ,M2 , VOL  2 , XP , XPROO , UE , QTOT , ( TSAMP! I >,1  =  1.9) 

190  FORMAT (9E14.6/9E14.6/9E14.6/) 

IF(NN.EQ.O)  GO  TO  87 
IPRINT  =  0 

90  1  PR  I NT= I  PR  I NT+ 1 
GO  TO  10 

99  CONTINUE 


STOP 

ENO 
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A  sample  printout  of  the  STG  program  follows.  Definitions  of  "Initial 
Condition"  variables  are  given  in  the  preceding  program  printout.  However,  some 
explanation  of  the  "Tabulated  Output"  variables  is  required. 


TIME 


PCH 


V0L1 


PI,  Tl,  Ml 

V0L2 


P2 ,  T2 ,  M2 
VELP,  XP 


Time  after  piston  release  -  seconds. 

TIME  is  initially  zero  and  increases  to  a  maximum  value 
either  when  the  projectile  exits  the  barrel  or  when  it 
equals  a  limiting  value,  i.e.,  TFO, 

Driver  gas  pressure  -  psia. 

PCH  is  a  prescribed  maximum  at  TIME  zero  and  decreases 
as  the  driver  gas  displaces  the  unlatched  piston. 

Test  gas  volume  included  from  piston  face  to  test 
specimen  inlet  -  ft3. 

V0L1  is  initially  the  entire  volume  of  the  driven  tube 
but  decreases  to  the  volume  of  the  test  gas  collection 
chamber  when  the  piston  has  been  fully  displaced. 

Test  gas  pressure,  temperature  and  mass  associated  with 
V0L1,  measured  in  psia,  °R  and  lbm,  respectively. 

Test  gas  volume  included  from  test  specimen  inlet  to 
projectile  base  -  ft3. 

V0L2  has  a  minimum  value  of  the  test  specimen  bore 
volume  at  TIME  zero  and  increases  with  projectile 
displacement  to  include  the  entire  barrel  volume. 

Test  gas  pressure,  temperature  and  mass  associated  with 
V0L2 ,  measured  in  psia,  °R  and  lbm,  respectively, 

Piston  velocity  and  displacement  -  ft/sec,  ft. 


VPROJ ,  XPROJ  Projectile  velocity  and  displacement  -  ft/sec,  ft. 

WORK  Work  performed  by  driver  gas  on  the  driven  piston,  and 

by  test  gas  on  the  projectile.  During  piston  rebound, 
the  program  also  computed  negative  work  done  on  the 
piston  by  the  test  gas  -  ft-lbf, 

UE  Test  gas  free  stream  velocity  -  ft/sec^. 

2 

QFLUX  Test  sample  surface  heat  flux  -  Btu/ft  -sec. 

2 

QTOT  Total  integrated  heat  input  to  test  sample  surface  -  Btu/ft  . 

TSAMP  Test  sample  surface  temperature  -  °R. 

TSAMP(l)  is  the  computed  temperature  on  the  test  sample  surface. 
TSAMP (2)  through  TSAMP (9)  are  subsurface  temperatures  computed 
at  depths  printed  in  the  nonlinear  DELX  array, which  is  amended 
when  the  time  increment  DELT  is  changed. 
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INPUT  OAT A 


TEST  NO. 


1 


TEST  CAS  MIXTURE  •¥  MOLE  TRACTION 

0.45S  Ml  THOSE* 

0.646  ARSON 


INITIAL  CONDITIONS  FOR  TEST 


P20  - 

2^  18.79980  PSF 

PCH  ■ 

91972.0000  PSF 

AT  • 

0.306800  FT«*2 

T20  ■ 

630.000000  0E6R 

TO  • 

636.000000  OECR 

KP1  • 

8.000000  L8F/FT 

R2  - 

46.173913  FT-tlf /L8M-P6CR 

RDVR  - 

66. 199997  FT-L8F/L0M-D66R 

USHELL  • 

160.000000  LBM 

CV2  - 

0.121737  BTU/LBM-OECR 

CVDVR  • 

0.176000  OTU/LftM-OEfiR 

SHELRV  - 

160.000000  L6M 

CVT2  • 

0.000669 

VOIQ  - 

31.600000  FT**3 

XPMAK  • 

81  .000000  FT 

CVEXP  • 

0.0 

P START  • 

7200.00000  PSF 

ALFA2  « 

697.913330 

8ARL  • 

16.000000  FT 

8ETA2  ■ 

0.017169 

A2  - 

0.306600 

8ORC0  - 

0.099400  FT 

CAM  2  • 

1.487621 

03  • 

0.600000 

DOR 12  - 

0.041667  FT 

C02  - 

0.760000 

RESO  - 

7200.00000  LBM 

DEL  T  * 

0.001000  SEC 

VOL  20  • 

0.0  FTmm3 

RESS  • 

4320.00000  LIM/fT 

PRC  I  • 

10.000000 

VOL2F  • 

0.002230  FT**3 

RESC  - 

7200.00000  LRM 

TFO  • 

0.260000  SEC 

VOL  If  • 

0.082000  FT**3 

WPROJ  - 

0.260000  IftM 

AFtAST 

06LX- 


6 1 0 1 62970E -03 .  .7932113951-03,  .1031 1 7409E -02 i  .  1 34052646E -02 .  . 1 74260196E -02 .  . 226649600E -02 .  . 2945 1 2960C -02 , 


.  39286647 9E -02, 
.2402416241-01  . 
.160747180 
.946911706 
6 . 93542385 

6EN0 


4977269711-02. 
. 3 1 23 1 3922E-0 1 , 
.196971191 
1.22969311 
7.71604633 


.647043064E-02. 
. 40600794  IE-01 , 
.254762352 
1.S9859799 
f 0 . 0308642 


.8411662761-02, 
.627909861C-0I  , 
.331190826 
2.07816316 
13.0401020 


. 109360086E-OI. 
. 606 1 S 1 99 1 C-0 1 
.430647774 
2.70180861 
16.9621179 


.1421649921-01. 
.  8919966226-01 . 
.669711764 
3.61208973 


,KTtMP» 


.  18480 1 370E  - 
.116969466 
.727624034 
4 . 6bb7 I  381 


01  , 


12. TJJP-  516.000000 


TABULATED  OUTPUT 


TIME 

PCH 

TSAMPI 1 ) 

PI 

P2 

TSAMPI 2) 

T 1 

T2 

TSAMPI 3) 

Ml 

M2 

TSAMPI 4) 

,  VOL  1 

1  VOL  2 

TS/yiP<6> 

VELP 

XP 

TSAMP 1 6  > 

VPROJ 

XPROJ 
TSAMPI 7) 

WORK 

uc 

TSAMPI 91 

OF  LUX 

OTOl 

TSAMPI  III 

0.0 

0 . 6352 1 3E  *03 
0. 8350006*03 

0 . 1 47774E*02 
0.1477746*02 
0.6350006*03 

0 . S306S4E  *03 
0.6306646*03 
0.6350006*03 

0.2166646*01 

0.1929126-03* 

0.5350006*03 

0.2493016*02 

0.2230006-02 

0 . 535000E*03 

0.5906746*01 

0.9099116*02 

0.6350006*03 

0.0 

0.0 

0.6360006*03 

0.0 

0.0 

0.5350006*03 

0.0 

0.0 

0.635000E*03 

0. 1 10000E-01 

0 . 632006E  *03 

0 .6350006*03 

0.1409716*02 

0.1499716*02 

0 . 636000E  *03 

0.6310946*03 

0.5310946*03 

0.6350006*03 

0.2166636*01 

0.193989E-03 

0.5350006*03 

0.2480446*02 

0.2230006-02 

0.5360006*03 

0.6417676*02 

0.9056156*02 

0.6360006*03 

0.0 

0.0 

0.6350006*03 

0.0 

0.0 

0.6350006*03 

0.0 

0.0 

0.6350006*03 

0. 210000E-OI 

0 . 623933E  >03 

0 . 636000E  *03 

0.1616806*02 

0.1616006*02 

0.6360006*03 

0.5349626*03 

0.6349626*03 

0.6360006*03 

0.21 5662E  *01 
0* 1962796-03' 
0 . 535UOOE  *03 

0.2460226*02 

0.2230006-02 

0.6360006*03 

0.1211726*03 

0.7959666*02 

0.6360006*03 

0.0 

0.0 

0.6360006*03 

0.0 

0.0 

0.6360006*03 

0.0 

0.0 

0.5350006*03 

0.31000UE-01 

0 . 61 1 4 1 9C  *03 
0.1360001*03 

0. 16S999E*02 
0.1669996*02 

0 . 63S000E  *03 

0 . 639878E  *03 
0.6390706*03 
0.6360006*03 

0.2156606*01 
0. 2001506-03 
0.6350006*01 

0.2402826*02 

0 .2230006-02 
0.5360006*03 

0.1764636*03 

0.7806146*02 

0.6350006*03 

0.0 

0.0 

0.6350006*03 

0.0 

0.0 

0.6360006*03 

0.0 

0.0 

0.6360006*03 

0. 409999E-01 
0.5960726*03 
0. 635000E *03 

0 . 1 623 1 2E  *02 
0.1623126*02 
0.6360006*03 

0 . 646732E  *03 
0.6467326*03 
0.5350006*03 

0.2166596*01 

0.2056246-03 

0.5350006*03 

0.2338836*02 

0.2230006-02 

0.6360006*03 

0.2296316*03 
0. 769659E *02 
0.6360006*03 

0.0 

0.0 

0.6360006*03 

0.0 

0.0 

0.6360006*03 

0.0 

0.0 

0.5350006*03 

0 . 609999E-OI 
0.8796101*03 

0  636000C *03 

0. I70013E*02 
0.1700136*02 

0 . 636000E  *03 

0.5656666*03 

0.6666666*03 

0.6360006*03 

0.2156676*01 

0.2128916-03 

0.6360006*03 

0.2258976*02 

0.2230006-02 

0.6360006*03 

0.2803816*03 

0.7338296*02 

0.6350006*03 

0.0 

0.0 

0.6350006*03 

0.0 

0.0 

0.6360006*03 

0.0 

0.0 

0.5350006*01 

0.609999E-01 

0.6637996*03 

0.6360001*03 

0.1019296*02 

0.1019296*02 

0.6360006*03 

0.6660076*03 

0.6660076*03 

0.6360006*03 

0.2156556*01 

0.2222296-03 

0.6350006*03 

0.2164036*02 

0.2230006-02 

0.6360006*03 

0.3284616*03 

0.7026836*02 

0.6360006*03 

0.0 

0.0 

0.6360006*03 

0.0 

0.0 

0.5360006*03 

0.0 

0.0 

0.5350006*03 

0 . 709999C -01 
0.6303496*03 
0.6360006*03 

0.1962076*02 

0.1962076*02 

0.6360006*03 

0.6006936*03 

0.6906036*03 

0.6360006*03 

0.2166536*01 
0.2340336  03 
0.6360006*03 

0.2064066*02 

0.2230006-02 

0.5360006*03 

0.3736696*03 

0.8871016*02 

0.6360006*03 

0.0 

0.0 

0.6360006*03 

0.0 

0.0 

0.636000E*03 

0.0 

0.0 

0.6360006*03 

0.6069606*03 

0.0360006*03 

0.2147996*02 

0.2147996*02 

0.1360006*03 

0.6974636*03 

0.6974636*03 

0.6360006*03 

0.2156606*01 

0.24B0656-03 

0 . 6350006  *03 

0.1932376*02 

0.2230006-02 

0.8350006*03 

0.4 168826 *03 
0.8271786*02 
0.6360006*03 

0.0 

0.0 

0.6360006*03 

0.0 

0.0 

0.0 

0.0 

0 .6350306*03 

80 


0  .  8099974  *  0 1  0.2368214*02  0 . 6  l  7762E *03 

0 . 4812344*03  0.2311211*02  0.6177624*01 

0 .6350004 ♦ 03  0.1310001*03  0.«36000E*03 

0 .  I  0 1 000E  *00  0.2704241*02  0.6423614*03 

0 ■ 4564344*03  0.2704241*02  0.6423614*03 

0 . S36Q00E  *03  0.6360004*03  0.6360004*03 

0.1  1  1000E*00  0.31 2l99f  *02  0.4724O2E*O3 

0 . 4  32BS 1E*03  0.3120994*02  0.6724U2£*Q3 

0 .6160004 *03  0.6360004*03  0.I3B000C*0J 

0.121  OOQE  *00  0.371  7341*02  0.7096424*03 

0 . 409274C  *03  0.371  736E*02  0.709642£*03 

0 . S3SOOOE  *  03  0 . 434000C  *03  0.4340004*03 

0 .  1  304994  *00  0.444704C»02  0.744374E*03 

0 . 38  70 1 2E  *03  0.4567041*02  0.766374E*03 

0.5350004*03  0.5360004*03  0.6340004*03 


0 . 2 1 664 7E  *0 1  0.179749C*02  0.4450524*03 
0.2676344-03  0.2230004*02  0.6932124*02 
0.5350004*03  0.6360004*03  0.6360006*03 

0.2156444*01  0.1641  164*02  0.49II38E»03 
0.2912414-03  0.2230004-02  0.6366164*02 
0.6350001*03  0.6360004*03  0.6360004*03 

0 . 1 1 66391 *01  0.1494334*02  0.624123E*03 
0. 32I801C-03  0.223000402  0.494396E«02 
0 . 5360001*03  0. ,6360004*03  0.6360004*03 

0.2166344*01  0.1327964*02  0.6639664*03 
0.3621114-03  0.2230004-02  0.430166E*02 
«. 6360004*03  0.6360004*03  0.636000E*03 

0.2166274*01  0.1162994*02  .0.6606614*03 
0.41 70444-03  0.2230004-02  0.3731364*02 
0.6360004*03  0.6360004*03  0.636000E*03 


0.0 

0.0 

0.636000E«03 

0.0 

0.0 

0.6360004*03 

0.0 

u.o 

O.SlbQuGF*03 

0.0 

0.0 

0.6360004*03 

0.0 

0.0 

0.6360004*03 

0.0 

0.0 

0 . 5150G0F.  *03 

0.0 

0.0 

0.6360004*03 

0.0 

0.0 

0.1)60004*03 

0 . 0 

0.0 

0.6350004*01 

0.0 

0.0 

0.6350004*03 

0.0 

0.0 

0 . 61&0004  *0  3 

0.0 

0.0 

O.bJbQOOt  *01 

0.0 

0.0 

0.6360004*03 

0.0  • 

0.0 

0.5350004*03 

0.0 

0.0 

0.b35000F«n > 

0.1409994*00  0.6967604*02  0.1171624*03  0.2166164*01  0.9704604*01  0.603973E*03  0.2962464*01  -0.2460624-03  0.0 

0.3669004*03  0.6667601*02  0.6171624*03  0.6047414-03  0. 2271774-02  0.3136444*02  0. 6496664-02  0.0  O.o 

0.6360004*03  0.6360004*03  0.1360004*03  0.6360004*03  0.6310004*03  0.6360004*03  0.6360004*03  0.6360004*01  0.5360004*03 

0.1609994*00  0.8030764*02  0.9994194*03  0.2166774*01  0.7914694*01  0.6236074*01  0.2670414*02  -0.4163724-02  0.0 

0.3460294*03  0.8030714*02  0.6994994*03  0.9941294-03  0.3204944-02  0.2620404*02  0.1262694*00  0.0  0.0 

0.6360004*03  0.6360004*01  0.6360004*03  0.6360004*03  0.6360004*03  0.6360004*03  0.6360004*03  0.6360004*03  0.535030F«03 

0.1609994*00  0.1211764*01  0.1019434*04  0.2164114*01  0.6073094*01  0.6302604*03  0.7796304*02-0.2313664-01  0.0 

0.3274691*03  0.121  1764*03  0.1010434*04  0.263904C-02  0.6922604-02  0.1M760E*02  0.6173824*00  0.0  0.0 

0.6360004*03  0.6360001*03  0.6360004*03  0.6360004*03  0.636000(*01  0.1110004*03  0.6360004*03  0.6360004*03  0.5350004*03 

4SECONO 

DELX-  .1929504664-03,  .2601364374-03.  .3260861364-03,  .4219112604-03,  .6610842994-03,  .7164099474-03,  .9313309574-03. 

. 121072936E-02.  .  1673947034-02,  .2046129904-02.  .2669967174-02,  .3467966004-02.  .4496337614-02,  .  684393367E - 02 . 

.7697107444-02.  .9076232694-02,  .4060078414-01.  .6270099601-01  .  .6061619014-01  .  .0919966221-01  ,  .  115969466 

.160747100  ,  .196971191  ,  .254762362  ,  .331190021  .  .430647774  ,  .669711764  .  .727624034 

.945911706  .  1.22960300  ,  1.59869799  ,  2.07916316  .  2.70160961  ,  3.61200973  ,  4.66671390 

5.93642305  7.71604633  .  10.0309642  13.0401020  16.0621179  .KTEMP-  16.TZIP-  667.743164  ,A- 


.0  .6-  .0 

64  NO 

0.1655994*00  0.1624744*03  0.1112304*04 
0.3167204*03  0.1296064*03  0.9790034*03 
0.6462944*03  0.6419434*03  0.6190404*03 

0 . 1 66599E  *00  0.1726204*03  0.1132914*04 
0. 3150264*03  0.1339374*03  0.9«00!3E*03 
0.6478304*03  0.6440014*03  0.6406204*03 

0.1676994*00  0.1039004*03  0.I16472E*04 
0.3133464*03  0.1309704*03  0.0034604*03 
0.6606104*03  0.6463364*03  0.6423014*03 

0 . 1 606094 *00  0.1964904*03  0.1177934*04 
0.3116004*03  0.1446701*03  0.9090094*03 
0.6633734*03  0.1497314*03  0.6441034*03 

0.1691994*00  0.2106144*03  0.1202704*04 
0.3100204*03  0.1007724*0)  0.9006204*03 
0.6664674*03  0.6612074*03  0.6461074*03 

0.1706994*00  0.2266324*03  0.1229214*04 
0.3093914*03  0.1676104*03  0.1000004*04 
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0.2904444*03  0.4314944*04  0.4499694*04 
0.2794174*04  0.2542994*04  0.2494734*04 

0.1922774*00  0.1419924*04  0.4092714*04 
0.2907234*03  0.4304074*04  0.4590276*04 
0.2740434*04  0.2430994*04  0.24441SE*04 

0.1923774*00  0.1643736*05  0.3999654*04 
0.2907944*03  0.4291444*04  0.4669644*04 
0.2696994*04  0.2497904*04  0.2449236*04 

0.1924774*00  0.1407244*05  0.3900344*04 
0.2909704*03  0.4249374*04  0.4634106*04 
0.2644064*04  0.2464346*04  0.2427624*04 

0.1926774*00  0.1377234*06  0.3914934*04 
0.2909496*03  0.4206921*04  0.4507466*04 
0.2612296*04  0.2530796*04  0.2404956*04 

0.1926776*00  0.1261916*05  0.3733416*04 
0.2910296*03  0.4169376*04  0.4479086*04 
0.2671766*04  0.2497556*04  0.2391376*04 

0.1927776*00  0.1169246*06  0.39S573E*04 
0.291  1  126*03  0 . 4 1 0422E  *04  0.4449356*04 
0.2632616*04  0.2464976*04  0.2357504*04 

0.1929776*00  0.1067966*05  0.3591776*04 
0.2811996*03  0.4045446*04  0.441656E*04 
0.2494936*04  0.2432936*04  0.2333506*04 

0.1929776*00  0.9866594*04  0.3511366*04 
0.281 7864  *03  0.3982934*04  0 . 4 38696E *04 
0.2458724*04  0.2401134*04  0.2309586*04 

0.1930774*00  0.9137911*04  0.3444346*04 
0.2013774*03  0.3917504*04  0.4354764*04 
0 . 2*23974*04  0.2371654*04  0.2295886*04 

0.1931776*00  0.8485786*04  0.3380526*04 
0.2814696*03  0.3849814*04  0.4322126*04 
0.2390664*04  0.2342416*04  0.2262516*04 

0.1932776*00  0.789682E*04  0.3319366*04 
0.2815626*03  0.3789676*04  0.4289816*04 
0.2358946*04  0.2314216*04  0.2239586*04 

0.1933766*00  0.7365816*04  0.3261046*04 
0.2816586*03  0.3728286*04  0.4267236*04 
0.2328646*04  O.2207I1E*O4  0.2217216*04 

0.1934766*00  0.6986976*04  0.3205476*04 
0.2817556*03  0.3663916*04  0.4224386*04 
0.2299466*04  0.2260906*04  0.2195386*04 

0.1935766*00  0.64S3806*04  0.3152486*04 
0.2818536*03  0.3697046*04  0.4191416*04 
0. 2271256*04  0.223549E*04  0.2174046*04 

0.1936766*00  0.6060886*04  0.3101906*04 
0.2819626*03  0.3528646*04  0.4158406*04 
0.2243916*04  0.2210776*04  0.2153146*04 

0.1937766*00  0.57035ZE*04  0.3053576*04 
0.2920546*03  0.3468946*04  0.4126496*04 
0.2217366*04  0.2196716*04  0.2132676*04 

0 . 1 938766*00  0.5377656*04  0.3007366*04 
0 . 282 1 666*03  0.3389376*04  0.409271E*04 
0.2191556*04  0.2163246*04  0.2112696*04 

0 . 1 93976E *00  0.6079796*04  0.2963126*04 
0.2822696*03  0.3317606*04  0.4060146*04 
0 . 2 1 6640E  *04  0.2140316*04  0.2092866*04 

0.1940766*00  0.4806866*04  0.2920766*04 
0.2823626*03  0.3246636*04  0.4027936*04 
0.2141866*04  0.2117896*04  0.2073496*04 

0.1941766*00  0.4666206*04  0.2990126*04 
0.2924896*03  0.3176736*04  0.3996836*04 
0.2117986*04  0.2095936*04  0.2064446*04 

0.1942766*00  0.4326606*04  0.2941136*04 
0.2928926*03  0.3038406*04  0.3932966*04 
0.7071136*04  0.2063146*04  0.2017166*04 

0.1944796*00  0.3916041*04  0.2767706*04 
0.2927906*03  0.2996166*04  0.3901936*04 
0.2049206*04  0.2032176*04  0.1999966*04 

0.1946766*00  0.3733946*04  0.273309E*04 
0.2929996*03  0.2997636*04  0.3971416*04 
0.7026326*04  0.2011336*04  0.1980896*04 


O  t^a  1  /6*01  0.1803046*00  -0.2688506*03 
0.1981866*00  0.7230656-01  0.2226336*00 
0.2231706*04  0.1984066*04  0.1498326*04 

0.1954126*01  0 . 199960E *00  -0.2907266*03 
0.2022316*00  0.7337246-01  0.2508486*00 
0.2238746*04  0.1908086*04  0.1526416*04 

0.1950416*01  0.1692496*00  -0.3105366*03 
0.2059436*00  0.7447846-01  0.2911256*00 
0.2241136*04  0.1922046*04  0.1650076*04 

0. 1947006*01  0.1711116*00  -0.3294886*03 
0.2093506*00  0.7862426-01  0.3132706*00 
0 . 223959E  *04  0.193526E*04  0.1572266*04 

0. I94396E*0I  0.1994934*00  -0.3447776*03 
0.2124616*00  0.7680926-01  0.3471096*00 
0.2234776*04  0.1946116*04  0.1592036*04 

0.1940976*01  0.1993476*00  -0.3596956*03 
0. 2153646*00  0.7I0333E-01  0.3924666*00 
0.2227306*04  0.1951966*04  0.1609446*04 

0.1939316*01  0.2108306*00  -0.3730906*03 
0.2190246*00  0.7929556-01  0.4192646*00 
0.2217676*04  0.1956176*04  0.I62462E*04 

0 . 1935846*01  0.222306E*00  -0 . 395407E *03 
0.2204856*00  0.9059536-01  0.4673216*00 
0.2206346*04  0.1959106*04  0.163769E*04 

0.1933566*01  0.234340E*00  -0 . 396699E *03 
0 . 222767E  *00  0.8193196-01  0.4965476*00 
0.2193686*04  0.1958064*04  0.I64882E*04 

0.1931434*01  0.2467046*00  -0 . 407069E *03 
0.2248886*00  0.0330426-01  O.S36846E*00 
0.2180006*04  0  .  1 95633E *04  0.1658166*04 

0.1929454*01  0.2593714*00  -0.4166194*03 
0.2260654*00  0.8471174-01  0.5781324*00 
0.2165576*04  0.1953164*04  0.1665864*04 

0 . 1 92760E ♦ 0 1  0 . 272317E  *00  -0.4254364*03 
0.2287114*00  0.8615334-01  0.6203284*00 
0.2I5060E*04  0.1948784*04  0.1672086*04 

0.1925874*01  0.285621E*00  -0 . 433597E *03 
0.2304404*00  0.0762014-01  0.6633664*00 
0.2135266*04  0.1943396*04  0.1676966*04 

0.1923746*01  0.2989646*00  -0 . 44  I  1 70E * 03 
0.2325654*00  0.891350E-01  0.707184E*00 
0.2119704 *04  0.I93714E*04  0.1680634*04 

0.1921686*01  0 . 3 1 2630E  *00  -0.44821  1E*03 
0.2346174*00  0.906733E-01  0.7517286*00 
0.2104054*04  0.193020E*04  0.I68323E*04 

0. 1 9 1 982E  *0 1  ■  0.3265036*00  -0.4547724*03 
0 . 2364  7  7  E  *00  0. 9224196-01  0.7969466*00 
0.2088416*04  0.1922704*04  0.1684866*04 

0.1910134*01  0.3405694*00  -0.4609004*03 
0.2301626*00  0. 9383994-01  0.8427944*00 
0 . 2072844  *04  0.1914754*04  0.1685654*04 

0. 19I660E*0I  0.3548164*00  -0.4666356*03 
0.2396054*00  0.9546664-01  0.8892314*00 
0.2057344*04  0.I90842E«04  0.1685676*04 

0. 191S226*01  0 . 369232E ♦ 00  -0.4720154*03 
0.2410584*00  0.9712076-01  0.9362194*00 
0 . 204 1 94E*04  0.1897784*04  0.1685024*04 

0.191 3984 *0 1  0 . 383B07E  *00  -0.4770704*03 
0.2422944*00  0.988016E-0I  0.9037274*00 
0.2026646*04  0.1888894*04  0.1683776*04 

0.1912074*01  0.3985316*00  -0.4818296*03’ 
0.2434006*00  0.1005086*00  0.103172E*01 
0.2011456*04  0.1079776*04  0.1681996*04 

0.1911896*01  0 .4133956*00  -0.4863176*03 
0.2443866*00  0.1022396*00  0.1080I7E«01 
0.1996376*04  0.1870486*04  0.167971E*04 

0.1911004*01  0.4283936*00  -0.4905576*03 
0 . 245259E  *00  0.103994E*00  0.1  129056*01 
0.1981416*04  0.1861026*84  0.1677016*04 

0.1910246*01  0.4435166*00  -0 . 494569E *03 
0.2466776*00  0.1075716*00  0.1226036*01 
0.1951806*04  0.1841736*04  0.1670446*04 

0 . 1 90902E *0 1  0.4741146*00  -0.5019776*03 
0.2472306*00  0.1093926*00  0.1278086*01 
0. 1937146*04'  0.1031926*04  0.1666656*04 

0.1908576*01  0.4895776*00  -0,5054046*03 
0.2476786*00  0.1112326*00  0.1329486*01 
0.1922836*04  0.1022016*04  0.1662576*04 


0.1376154*04  -0 . 1 39862E ♦ 00  0.1250834*05 
0.9219614*01  0.3162464*04  0.8053716*02 
0.1  134206*04  0.8548596*03  0.6*30404*03 

0.1426954*04  -0.1440556*00  0.1117584*05 
0.9359866*01  0.3116886*04  0.8171336*02 
0.1187616*04  0.8687886*03  0.6097464*03 

0.1481356*04  -0.1477146*00  0.9991176*04 
0.9505366*01  0.3071926*04  0-82/6466*02 
0.1190196*04  0.8829816*03  0.6959156*03 

0. 1533446*04  -0.1508594*00  O.S94464F*04 
0.9656104*01  0.3028014*04  0.8370486*07 
0.1202094*04  0.8973526*03  0.7023444*03 

0 . 1  58494E  *04  -0 . 1 535 1 SE *00  0.0023556*04 
0.9812026*01  0 . 290542E  *04  0.8454706*02 
0.1223064*04  0.9118176*03  0.7090214*03 

0 . 1 63560E *04  -0.1557086*00  0.72I429F*04 
0.9973064*01  0.2944316*04  0.8530456*02 
0 . I 24  302E  *04  0.9262946*03  0.7159304*03 

0.1685566*04  -0.1574706*00  0.650353E*04 
0.1013916*02  0.2904776*04  0 . &V.98C.0F  *0? 
0.1261884*04  0.940/046*03  0  . /230b24  «  03 

0 . 1  734 45E  *04  -0.1598294*00  0.5878864*04 
0 . 1 03 1 01 E*02  0.2866824*04  0.8660124*02 
0.1279606*04  0.9549756*03  0.7303666*03 

0.1782274*04  -0.1598146*00  0.53251*'!:: 
0.1040606*02  0.2830466*04  0.8715824*02 
0.1296176*04  0.9690404*03  0 . 7  .’OS  It  *  03 

0.1828964*04  -0 . I 60456E • 09  0.4044346*04 
0.1066656*02  0.2795654*04  0 . 8 /6b39E  *02 
0.131  1584*04  0 . 90284 1 E  *03  0.74541121*03 

0.1874454*04  -0.1607826*00  0.44IS*5E*04 
0.1085174*02  0.2762336*04  0.0812424*02 
0.1325844*04  0 . 9O6320E ♦ 03  0.7632354  *03 

0.1918704*04  -0.  1608214*00  0.4036?»4*04 
0.1104144*02  0.2730456*04  0.88^4454*02 
0 . 1  33899E  *04  0.1009464*04  0.76108*4*03 

0.1961696*04  -0.1605984*00  0.3698344*0/ 
0.1123546*02  0.269993E*04  0 . 089292f ♦ 02 
0.1351074*04  0.1022204*04  0 . 76900  f • 03 

0.2003444*04  -0.1604956*00  0 . 34 1 46 ?F *04 
0 . I  1 4337E  *02  0.2584346*04  0  8928294*02 
0.1362114*04  0.1034524*04  0.7769786*03 

0.2044024*04  -0.1602334*00  0.3153/36*04 
0  .  1  1 6360E  *02  0.2442904*04  0.8960984*02 
0.1372166*04  0.1046404*04  0.7849746*03 

0.2083404*04  -0.1596906*00  0.2908224*04 
0.1184246*02  0.2304246*04  0.0991156*02 
0.1381296*04  0.1057826*04  0.792973F*03 

0.2121546*04  -0.1588986*00  0.26/7206*14 
0 . 1  2052  7E  *02  0.2167886*04  0  .  «>0 1  895F  *  >2 
0.1389536*04  0.1068784*04  0.8009556*03 

0.2158434*04  -0.1578846*00  0.2459586*04 
0.1226676*02  0.2033354*04  Q.y04/50f«02 
0.1396964*04  0.1079286*04  0. 0089096*03 

0.21 9406E *04  -0.1566736*00  0.2254096*04 
0.1248434*02  0.1900066*04  0 . 906 7y.t  ♦  02 
0.1403614*04  0.1089306*04  0.8167916*03 

0.2228434*04  -0 . 1552906*00  0 . 20593',E *04 
0 . 1270546*02  0.1767386*04  0.9089414*02 
0  .  I  40954E  *04  0.1096864*04  0 . 82460°E *03 

0.2261556*04  -0.1537566*00  0.187412E*04 
0.1292996*02  0.1634586*04  0.910897F*02 
0. 1414816*04  0.1  107964*04  0.8323424*03 

0 .2293416*04  -0.1520926*00  0.1696734*04 
0.1315776*02  0.1500746*04  0.9126776*02 
0.1419446*04  A.!11660E*04  0.8399746*03 

0.2324044*04  -0.1503136*00  0.1525526*04 
0.1338064*02  0.1364676*04  0.9142736*02 
0.1423494*04  0.1124816*04  0 . 84 7495F *03 

0 . 23S346E  *04  -0.1484386*00  0.1358304*04 
0.1385926*02  0.1078666*04  0.9169726*07 
0.1429996*04  0.1139036*04  0.8621594*03 

0.2408706*04  -0.1444484*00  0. 1023704*04 
0.1409876*02  0.9224704*03  0.9100/1692 
0.1432524*04  0.1  146896*04  0  .  *0.1 

0.2434674*04  -0 . 1 42355E *00  0.8441886*03 
0.1434096*02  0.7490414*03  0.91*9976*02 
0.1434596*04  0.1163434*04  0.0762694*03 
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APPENDIX  II 


ISENTROPIC  EQUILIBRIUM  COMBUSTION  CODE 


Adiabatic  Compression  Program 

An  existing  equilibrium  combustion  program  has  been  modified  to 
calculate  the  properties  of  a  gas  mixture  at  various  points  on  an  isentrope. 
The  program  uses  a  list  of  thermal  properties  of  possible  constituents  of  the 
mixture.  The  program  writes  output  containing  the  composition  and  thermal 
properties  of  the  gas  mixture  at  the  input  pressures. 

Input  to  the  program  consists  of  "REACTANT"  cards  specifying  the 
composition  at  the  initial  point  on  the  isentrope  and  a  namelist,  "SINPT2," 
which  specifies  the  temperature  and  pressure  at  the  initial  point  and  a 
schedule  of  other  pressures  (up  to  26  pressures  including  the  initial  value) . 
OMIT  cards  may  be  used  if  desired  to  exclude  certain  species  from  considera¬ 
tion  as  mixture  constituents. 


The  program  operates  as  follows: 


(1)  The  initial  mixture  as  specified  by  the  reactant  cards  and  the 

mixture  ratio  specified  by  §INPT2  is  analyzed  to  determine  which 
chemical  elements  are  present.  The  number  of  gram-atoms  of  each 
element  is  calculated  as  well  as  the  mixture  molecular  weight 
and  mixture  total  enthalpy  (Subroutine  REAC) . 


(2)  The  taped  thermal  data  is  searched  (SEARCH)  and  the  names  and 

thermal  properties  of  possible  compounds  which  could  be  formed 
from  the  available  atoms  are  extracted.  This  list  is  then  com¬ 
pared  to  the  OMIT  cards  and  the  net  list  of  species  to  be 
considered  is  printed. 


(3)  Subroutine  EQLBRM  then  varies  the  composition  seeking  to  minimize 

the  Gibbs  free  energy  of  the  mixture  subject  to  the  constraint 
that  the  numbers  of  atoms  of  the  different  elements  do  not  vary. 


In  this  process,  species  with  molefractions  of  less  than  10  at  any 
iteration  are  dropped  from  consideration.  Occasionally,  the  constraint 
equations 


where  a^  is  the  number  of  atoms  of  type  i  in  species  j,  become  linearly 
dependent.  In  this  case,  the  program  will  print  the  message  "SINGULAR  MATRIX" 
and  then  calculate  the  thermal  properties  of  the  input  mixture  at  the  specified 
pressure  and  temperature  (Subroutine  HCALC) .  Sometimes  the  mixture  free  energy 
will  have  a  very  broad  minimum  or  will  possess  two  minima 
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as  in  the  sketch.  If  this  is  the  situation,  the  program  will  print  the 
message 


"35  ITERATIONS  DID  NOT  SATISFY  CONVERGENCE  REQUIREMENTS  FOR  THE 
POINT  _ " 

and  then  calculate  the  thermal  properties  of  the  input  mixture. 

(4)  One  way  or  the  other,  the  program  will  obtain  the  entropy  at  point  #1, 
This  entropy  will  then  be  held  constant  (SP=.TRUE.).  EQLBRM  will 
calculate  the  equilibrium  composition  at  an  estimated  temperature 

and  the  new  pressure,  calculate  the  mixture  entropy,  and  compare  it 
to  the  entropy  at  point  #1.  If  the  two  entropies  do  not  agree,  the 
temperature  estimate  is  changed  and  a  new  equilibrium  calculated. 

(5)  If,  at  any  point,  the  equilibrium  calculation  fails,  the  control 
program  THERMP  will  call  subroutine  FR0ZEN  to  obtain  the  thermal 
properties  corresponding  to  the  fixed  entropy,  the  new  pressure 
and  the  composition  at  the  previous  point. 

(6)  After  13  points  have  been  calculated,  the  program  will  print  the 
results.  Then  it  will  complete  and  output  the  remainder  of  the 
pressure  schedule, 

(7)  After  completing  the  pressure  schedule,  the  program  will  look  for 
a  new  mixture  ratio  MIX. 


REACT.*  NTS 

AR  l.y.w 

0.0 

0.0 

0.0  oo 

0.0 

0.1410 

M 

0.0 

c 

290.100  F 

0.0 

N  2.0000 

O.C 

0.0 

0.0  00 

0.0 

0 . 4560 

M 

0.0 

G 

291. 150  F 

0.0 

OMIT 

NAMELISTS 

CIS) 

AINPT2 

•CASE- 

2 ,  T- 

.  .0 

.  .0 

.  .0 

.  .0 

.  .0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.6 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

-0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

14. 

26.0000000 

40.0000000 

, 

134.000000 

446.000000 

1340.00000 

2233. 

)0000 

4400. 

10000 

•699.00000 

•932.00000 

11165.0000 

13391.0000 

15631.0000 

17104 

.0000 

20097 

.0000 

22330.0000 

21603.0000 

26796.0000 

29029.0000 

31262.0000 

I340I 

.0000 

30721 

.0000 

37961.0000 

40104.0000 

42427 . 0060 

44660.0000 

,  PS  1 A»T , MMHG-F , NSOM-F , V- 

.0 

, 

.0 

.  .0 

.  .0 

, 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.  0 

.0 

,  21.0000000 

,  41.0000000 

,  134.000000 

,  446 . OOOOOC 

1340 .00000 

2233.00000 

4466.00000 

0099 . 

10000 

0932. 

>0000 

1 1166.0000 

13300.0000 

15631.0000 

17164.0000 

20097.0000 

22130 

.0000 

24103 

.0000 

20796.0000 

29029.0000 

31262.0000 

33495.0000 

36720.0000 

37961 

.0000 

40194.0000 

42427.0000 

44000.0000 

.ERATIO-F 

.  .FPCT-f ,FA-F 

MIX- 

.0 

.  .0 

.  .0 

,  .0 

.0 

.  .0 

.0 

. 

.0 

.0 

.0 

♦ 

.0 

. 

.0 

.0 

* 

so- 

AEND 


, IONS"F , 1DIIUG" 


0, TRACI-  .0 


.IIUNIT-F .EUNITS 


F ,TRNSPT*F . FROZN-F .PUNCK-F .NODATA-F , 0 1 F -f 


NO  INPT2  /AlUE  GIVEN  FOR  OF,  IORAT,  FA,  OR  FPCT 
SPECIES  IE  I NC  CONSIDERED  IN  THIS  SYSTEM 


L  5/66 

AR 

J  3/77  N 

J  3/77  N2 

J 12/70  N3 

OF 

.  0  0 

EFFECTIVE  FUEL 

EFFECTIVE  ONIDANT 

MIXTURE 

ENTHALPY 

HPP( 2  > 

HPP(  1  1 

HSUIO 

(KG 

-  MOL  HO 

Kl/KG 

0.0 

0.0 

0.0 

KG- 

ATOMS /K  , 

OOP! 1,2) 

I0P( 1 , 1  1 

•0<  I  > 

AR 

0 . I 67S0972C-0 1 

0.0 

0.1670M72E-0! 

N 

0 . 26363246E -01 

0.0 

0 . 26363246E-01 

PT 

AR 

N 

1 

-19.217 

-11.910 

11.000 

2 

-10. OOO 

-11.646 

1.000 

3 

-10.223 

-11.430 

1.000 

4 

-17.56! 

-11.190 

2.000 

5 

-10.965 

-11.023 

2.000 

6 

-10.405 

-10.942 

2.000 

7 

-10.205 

-10.924 

2.000 

1 

-10.025 

-lO.fli 

2.000 

9 

-15.077 

-10.910 

1.000 

10 

-10.770 

-10.923 

1.000 

II 

-10.099 

-10.929 

1.000 

12 

-15.030 

-10.930 

1.000 

13 

-15.604 

-10.942 

1.000 

THERMODYNAMIC  EQUILIBRIUM  PROPERTIES  AT  ASSI6NE0 
ENTROPY  AND  PRESSURE 


CASE  NO.  2 

CHEMICAL  FORMULA 
FUEL  AR  1.00000 

FUEL  N  2.00000 


MOLES 

0.54500 

0.46500 


ENERGY  STATE 
CAL /MOL 

0.0  G 

0.0  G 


TEMP  DENSITY 
DEG  K  C/CC 
291.16  0.0 

291.15  0.0 


O/F-  0.0  PERCENT  FUEL-100.0000  EQUIVALENCE  RATIO- 


REACTANT  DENSITY-  0.0 


THERMODYNAMIC  PROPERTIES 


P.  ATM 

T.  DEG  K 

RHO,  G/CC 

H.  CAl/G 

S.  CAL  /( 6HKI 

1.0003 

300 

1.4020-3 

0.3 

1.2279 

1.7011 

369 

1 .9932-3 
10.4 

1 .2279 

3.0621 

430 

2.9414-3 

23.9 

1.2290 

9.1161 

631 

6.0027-3 

67.0 

1.2279 

30.340 

933 

1.3091-2 

110.0 

1.2270 

91.101 

1310 

2.0106-2 
lit  .1 
t . 2276 

161.00 

1842 

4.1466-2 

223.1 

1.2279 

303.09 

1902 

6.7204-2 

291.0 

1.2279 

456.94 

2149 

9.9227-2 

339.0 

1.2260 

607.76 

2341 

1.0921-1 

376.0 

1.2200 

759.73 

2601 

1.2777-1 

407.1 

1 .2279 

911.00 

2640 

1.4626-1 

434.0 

1.2279 

1063.62 

2764 

1 .6190-1 
450.0 

1 .2279 

M.  MOL  WT 

(DLV/0LP1T 

(DLV/0LT1P 

CP.  CAL/C  G  H  K 1 
GAMMA  (SI 

SON  VEL.H/SEC 

34.519 

-1.00000 

1 .0000 
0.1702 
1.0111 
330.4 

34.011 

-1.00000 

1.0000 

0.1703 

1  .6100 
301.4 

34.510 

-1.00000 

1.0000 

0.1709 

1.6000 

390.0 

34.611 

-1.00000 

1.0000 

0.1739 

1.4949 

476.6 

34.611 

-1.00000 

1.0000 

0.1002 

1.4604 

674.6 

34.119 
-1.00000 
1.0000 
0.1061 
1.4479 
•  76.0 

34.810 

-1.00000 

1.0000 

0.1006 

1.4396 

731.2 

34.610 

-1.00000 

1.0000 

0.1913 

1.4305 

•09.6 

34.619 

-1.00000 

1 .0000 
0.1926 
1.4263 
160.3 

34.610 

-1.00000 

I .0000 
0.1934 
1.4239 
•96.0 

34.810 

-1.00000 

1.0000 

0.1939 

1.4223 

926.7 

34.610 
-1 .00000 

1 .0000 
0.1943 
1.4211 
960.7 

34.610 

-1.00000 
1.0000 
0. 1946 
1.4202 
972.3 

MOLE  FRACTIONS 

.  • 

AR 

N2 

0.64600 

0.46600 

0.64900 

0.46600 

0.64600 

0.46600 

0.64600 

0.46600 

0.64600 

0.46000 

0.64600 

0.40600 

0.64600 

0.40600 

0.64600 

0.46600 

0.64600 

0.46000 

0.64000 

0.48600 

0.64800 

0.46600 

0.54600 

0.46600 

0.64600 

0.45500 

ADDITIONAL  PRODUCTS  WHICH  WERE  CONSIDERED  IUT  WHOSE  MOLE  FRACTIONS  WERE  LESS  THAN  0.50000E-08  FOR  ALL  ASSIGNED  CONDITIONS 
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r 


N 


l 

-  IB. SIS 

-10.940 

7.000 

2 

-  IB .  490 

-10.991 

1.000 

3 

-  IS. 403 

- 10. 90S 

1.000 

4 

-11.412 

-10.903 

1.000 

s 

-  IS . 403 

-10.900 

I.  000 

0 

-19.170 

-10.973 

1.000 

7 

-10.102 

-10.977 

1.000 

0 

-IS. 329 

-10.901 

2.000 

• 

- 10. 300 

—10. 900 

1.000 

10 

-10. 200 

-10.909 

i  .000 

11 

- 10 . 209 

-10.993 

1 .000 

12 

-19.202 

-10. 990 

1 .000 

13 

- 10.230 

-11.000 

1.000 

THERMODYNAMIC  EOUILI9RIUM  MOMITIH  AT  ASSIGNED 
ENTROPY  AND  MIIIUM 


NOUS  ENERGY  STATE  TINT  DENSITY 

CAL /MOL  DEC  K  fi/CC 

0.94900  0.0  6  ISO. IS  0.0 

0.4SS00  0.0  6  IDS. IS  0.0 

PERCENT  FUEL* 100 . 0000  EQUIVALENCE  RATIO-  0.0  REACTANT  DENSITY-  0.0 


CASE  NO.  2 

CHEMICAL  FORMULA 
rUEL  AR  1.00000 

FUEL  N  2.00000 

0/F-  0.0 


THERMODYNAMIC  PROPERTIES 


P,  ATM 

1219.07 

1397.01 

1019.40 

1071.41 

1923.39 

1970.30 

2127.24 

2279.19 

2431.14 

2093. 09 

2736.03 

2999.97 

3039.92 

T.  OIG  K 

2070 

2977 

3070 

3100 

3240 

3317 

3390 

3400 

3029 

3999 

3690 

3709 

3769 

RHO,  C/CC 

1.7710-1 

1.9327-1 

2.0917-1 

2.2204-1 

2.3074-1 

2.0040-1 

2.9394-1 

2.7712-1 

2.9004-1 

3.0273-1 

3.1921-1 

3.2749-1 

3.3909-1 

H.  CAL/G 

479.0 

499. S 

017.9 

034.9 

000.9 

000.0 

090.3 

093.0 

009.9 

619.3 

931.3 

942.7 

603.7 

S,  CAL/ ( C X  K ) 

1.2279 

1.2279 

1.2279 

1.2279 

1.2279 

1.2270 

1.2279 

1.2279 

1.2279 

1.2279 

1.2279 

1.2279 

1.2279 

M.  MOL  WT 

34.010 

34.010 

34. SIS 

34.019 

34.019 

34.011 

34.019 

34.010 

34.919 

34.919 

34.919 

34.019 

34.019 

(DLV/DLP1T 

-1.00000 

-1.00000 

-I. 00000 

-1.00000 

-1.00000 

-1.00000 

-1.00001 

-1.00000 

•1.00000 

-1 .00000 

-1.00000 

•1.00000 

-1.00000 

<  DLV/DLT1P 

1.0000 

1.0000 

1.0000 

1.0000 

1  .0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0001 

1.0001 

CP.  CAL/IGMK) 

0.1940 

0.1900 

0.1902 

0. 1904 

0.1900 

0.1000 

0.190B 

0.1009 

0.1960 

0.1961 

0.1993 

0.1964 

0.1960 

GAMMA  <S> 

1.4194 

1.4100 

1 .4103 

1.4170 

1.4174 

1.4199 

1.4109 

1.4162 

1.4166 

1.4169 

1.4192 

1.4I4B 

1.4146 

SON  VCL .M/SEC 

991.4 

1000.9 

1024.2 

1030.0 

1051.7 

1064.0 

1070.0 

1080.4 

1096.6 

1106.2 

1119.4 

1124.2 

1132.5 

MOLE  FRACTIONS 


AR 

0.04800 

0.04000 

0.04000 

0.04500 

0.04000 

0.04500 

0.84600 

0.64500 

0.54000 

0.64600 

0.54500 

0.04000 

0.04600 

N 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00001 

0.00001 

0.00001 

N2 

0.45900 

0.48800 

0.46500 

0.49000 

0.45500 

0.45500 

0.46500 

0.45500 

0.45500 

0.40600 

0.40900 

0.40000 

0.45499 

ADDITIONAL  PRODUCTS  WHICH  WERE  CONSIDERED  BUT  WHOSE  MOLE  FRACTIONS  WIRE  LESS  THAN  O.SOOOOE-OS  FOR  ALL  ASSIGNED  CONDITIONS 
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Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet  and  return  it  to  Director,  US  Army  Ballistic  Research 
laboratory,  ARRADCOM,  ATTN:  DRDAR-TSB,  Aberdeen  Proving  Ground, 
Maryland  21005.  Your  comments  will  provide  us  with  information 
for  improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/ contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 


Telephone  Number: 
Organization  Address: 


